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PBEFACE. 



In March, 1875, the usual biennial notice was issued, giving the 
subjects for the Adams Prize to be adjudged in 1877. The 
following is the chief portion of the notice : 

The University having accepted a Fund raised by several members 
of St John's College for the purpose of founding a Prize to be called the 
Adami Frize, for the best essay on some suliject of Pure Mathematics, 
Astronomy or other branch of Natural Philosophy, the Prize to be 
given once in two years, and to be open to the competition of all 
persons who have at any time been admitted to a degree in thi^ 
University — 

The Examiners give notice that the following is the subject of the 
Prize to be adjudged in 1877 : The Criterion of Dymmicoil StabiMy. 

To illustrate the meaning of the question imagine a particle to slide 
down inside a smooth inclined cylinder along the lowest generating line, 
or to slide down outside along the highest generating line. In the 
former case a slight derangement of the motion would merely cause 
the particle to oscillate about the generating line, while in the latter 
case the particle would depart from the generating line altogether. 
The motion in the former case would be, in the sense of the question, 
stable, in the latter unstable. 

The criterion of the stability of the equilibrium of a system is, 
that its potential energy should be a minimum; what is desired is, a 
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corresponding condition enabling us to decide when a dynamically pos- 
sible motion of a system is such, that if slightly deranged the motion 
shall continue to be only slightly departed from. 

The essays must be sent in to the Yice-Chancellor on or before the 
16th December 1876, &c., <fec. 

S. G. PHEAK, Vtce-ChanceUor. 

J. CHALLIS. 

G. G. STOKES. 

J. CLERK MAXWELL. 



The pressure of other engagements for some time prevented 
me from giving my attention to the subject. This essay was 
therefore almost entirely composed during the year 1876. It is 
now printed as it was sent in to the Examiners, the changes being 
merely verbal. Some few additions have been made where ex- 
planation appeared to be necessary, but all these have been 
marked by square brackets, so that they can be at once dis- 
tinguished from the original parts of the essay. 

In order to shorten the essay as much as possible many 
merely algehraic processes have been omitted and the results only 
are stated. It is hoped that this will add clearness as well as 
brevity to the reasoning, as the attention of the reader will not 
be called from the argument to follow a manipulation of symbols 
which may not present any novelty. 

The line of argument taken may be indicated in a general way 
as follows. Chapter I. begins with some definitions of the terms 
stahle and steady motions. It is then pointed out that whether 
the forces which act on the system admit of a force-function or not, 
the stability of the motion, if steady, is indicated by the nature 
of the roots of a certain determinantal equation. The boundary 
between stability and instability being generally indicated by the 
presence ^f equar roots, a criterion is investigated to determine 
beforehand whetheT^qual roots do or do not imply instability. 
This case being disposed of, the consideration of the determinantal 
equation is resumed. Two general methods are given by which, 
without solving the equation, it may be ascertained whether the 
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character of the roots imply stability or instability. These occupy 
Chapters !!• and III. In the first method a derived equation is 
made use of, and it is shown that a simple inspection of the signs 
of the coeflScients of the several powers in these two equations will 
decide the question of stability. In the second method a certain 
easy process is found which if performed on the determinantal 
equation will lead to the criteria of stability. At the end of the 
third Chapter a geometrical interpretation is given to the argu- 
ment. 

In the fourth Chapter the forces which act on the system 
are supposed to have a force-function. The determinantal equa- 
tion is then much simplified. Several points are considered in 
this Chapter which are necessary to the argument, such as the 
proper method of choosing the steady co-ordinates (if there be 
any), the distinction between harmonic oscillation about steady 
motion and that about equilibrium, and the changes which must 
be made in the determinantal equation when the equations of 
Lagrange become inapplicable. A method of modifying the 
Lagrangian function is also given by which, in certain cases, 
the fundamental determinant may be reduced to one of fewer 
rows and columns. 

In the fifth Chapter a series of subsidiary determinants is 
formed, and it is shown that at least as many of the conditions 
of stability are satisfied as there are variations of signs lost in the 
series in passing from one given state to another. It is also shown 
that this is equivalent to a maximum condition of the Lagrangian 
function. 

In the sixth Chapter the energy test of stability is considered. 
It is also shown that, when the motion is steady, this reduces 
to the same criterion as that indicated in Chap. V. 

In the seventh Chapter the question considered is whether 
the stability of a state of motion can really be determined by an 
examination of the terms of the first order only. In some cases 
these are certainly sufficient, and an attempt is made to discrimi- 
nate between these cases and those in which the terms of the 
higher order ultimately alter the character of the motion. 

If the Hamiltonian characteristic and principal functions be 
given, the conditions of stability as regards space only, or both 
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space and time may be deduced. But if the^e be not known as 
expressed in the Hamiltonian form, we may yet sometimes dis- 
tinguish between stability and instability if we can determine 
whether a certain integral ceases to be a minimum at some instant 
of the motion. This is the subject of the eighth chapter. 

As part of the third edition of my treatise on the Dynamics of 
Bigid Bodies was written at the same time as this essay, there are 
necessarily points of contact between the two works. Thus the 
subjects of the first part of the seventh chapter and of a portion 
of the sixth will be found discussed in the treatise on Dynamics. 
But as the objects of the two books are not the same, it will be 
found that in all these ca^es there are considerable differences in 
the modes of demonstration* 

EDWARD J. ROUTH. 

Petbrhoube, 

August 14, 1877. 
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CHAPTEE I. 



Dejinitiona of the terms small quantity, stable motion, steady motion. 
Arts. 1, 2. 

A system of bodies in steady motion is stable if the roots of a certain 
determina/ntaZ equation a/re such that their real pa/rts a/re oM negative. 
Art. 3. 

Effect of eqtml roots, and a test to determine whetlier equal roots do or do 
not introduce terms which contain the time as a factor. Arts. 4 — 7. 

Object of Chapters II, and III. Art. 8. 

1. Let us suppose a dynamical system to be set in motion 
under any forces and to move in some known manner. If any 
small disturbance be given to the system, it may deviate only 
slightly from its known motion, or it may diverge further and 
further from it. Let 6, <f>, &c. be the independent variables or 
co-ordinates which determine the position of the system, and let 
the known motion be given by ^ = ^q, = ^^, &c. where d^y (fy^y &c. 
are known functions of the time t. To discover the disturbance of 
the system we put = 0^-\-x, = ^o + y» ^^- These quantities 
Xy y, &c. are in the first instance very small because the disturb- 
ance is small. The quantities a?, y, z, &c. are said to be small 
when it is possible to choose some quantity numerically greater 
than all of them, which is such that its square can be neglected. 
This quantity may be called the standard of reference for small 
quantities. 

If, after the disturbance, the co-ordinates x, y, «, &c. remain 
always small, the undisturbed motion is said to be stable; if, on the 
other hand, any one of the co-ordinates become large, the motion is 
called unstable. 

It is clear that the same motion may be stable for one kind of 
disturbance and unstable for another. But it is usual to suppose 
the disturbance general, so that if the motion can be made un- 
stable by any kind of disturbance (provided it be small) it is said 
to be unstable. On the other hand, it will be called stable only 
when it is stable for all kinds of small disturbances. 

K. A. 1 
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2 DEFINITIONS. [CHAP. 

2. To determine whether x, y, z, &c. remain small, we must 
substitute for 0, <]>, &c. in the equations of motion their values 
^0 ■♦■ ^> ^0 + y> ^^' Assuming that x, y, &c. remain small, we may 
neglect their squares, and thus the resulting equations will be 

linear in x, y, jg, &c. The coefficients of a?, -^ , -^ >y»~Jfy ^' *'^- 

in these equations may be either constants or functions of the 
time. In the former case the undisturbed motion is said to be 
stead y foTlKese' co-ordinates, m the latter wnstewdy. In the case 
ot a steady motion a?, y, «, «&c. are all functions of the time 
which has elapsed since the disturbance and of certain constants 
of integration which are determined by the initial values of 

dx dti 
a?, -^ , y, -^ , &c. We may therefore define a steady motion to be 

such that the same change of motion follows from the same initial 
disturbance at whatever instant the disturbance is communicated 
to the system. 

If all the coefficients in the equations to find x, y, z are con- 
stant, they may be made to contain ^ by a change of co-ordinates. 
Thus we may write for x, y, z. &c. 

y = a'f + /3'i7+... 

where a, /8, &c. are any functions of t we please. Conversely, when 
the coefficients are functions of t, we may sometimes make the 
coefficients constant by a proper change of co-ordinates. But this 
cannot always be done. If there are n co-ordinates, we have v? 
arbitrary functions a, )8, &c. at our disposal. In each of the n 
linear equations of motion we may have three terms for each 
co-ordinate, and thus we have (3 n — 1) w coefficients to make con- 
stants. We have therefore in general too many equations to 
satisfy. The proper method of choosing the co-ordinates of refer- 
ence will be considered in a future chapter. 

3. Let us suppose a dynamical system to be making small 
oscillations under the action of any forces which may, or may not, 
possess a force function and to be subject to any resistances which 
vary as the velocities of the parts resisted. The general equations 
of motion will then be of the form 

&c. = 
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1.] THE DETERMINANTAL EQUATION. 3 

To solve these equations we write 

aj = lfo"^, y = Jf'e'»*, &c. 

Substituting we obtain a determinantal equation to find m. 
If we put 

A^^A^rnJ' + A^m + A,, B^Bji^ + B^m^-B^, 
A! = ^>' + Aim + ^;, &c. 

this equation may be written in the simple form 



A, B, (7, 
A, B, G\ 



= 0. 



We may also write the equation in the form/(m) = 0. 

The coefficients M, M\ &c. are not independent, but if we 
represent the minors of -4, jB, (7, &c. by a, J, c, &c. we may easily 
show that 



We also have 



a b c 



M M M' , 
a o c 



It may be shown by properties of determinants that these 
equations all give the same ratios. If A, be the second minor 
obtained from the determinant f{m) by omitting the first and 
second rows and columns, we know that 

A,/ (m) = ah' - cih. 
Hence iff{m) = we have - = p • 

ft p 

In the same way we may show that -7 = - , and so on. This 

CL C 

property of Determinants is given in Dr Salmon's Higher Algebra, 
Lesson iv. Ex. 1. 

The general solution of the equation may therefore be written 
in the form 

X = Lfi^^^ + L^a^e^^ + . . . ] 

z = &c. J 

where ij, i,... are arbitrary constants, a^, a, &c. the values of 
the minor a when m^, m, ... are substituted for m ; ftp 6j ... the 

1—2 



4 NATURE OF THE STABILITY [CHAP. 

values of the minor b when similax substitutions are made, and 
so on. 

4. We see that the whole character of the motion will de- 
pend on the signs of the quantities m^, m,... If any one be real 
and positive, x, y, &c. or some of them will ultimately become 
large, and the steady motion about which the system is oscillating 
will be unstable. If all the roots are real, negative or zero and 
unequal, the motion will be stable. 

If two of the roots be imaginary we haVe a pair of imaginary 
exponentials. If these imaginary roots be a ± )3V — 1, the terms 
can be rationalized into 

e^(N^ cos I3t + N^ sin ^t). 

The motion will be stable if a be negative or zero, and unstable 
if a be positive. 

If two roots be equal, the form of the solution is changed. Let 
m, = m^ + A where h will be ultimately zero, we then have 

If we now make L. and Z, infinite in the usual manner, we 
find 

&c. = &c., 
where J^, M^ are two arbitrary constants which replace i^, L^. 

In the same way if three roots are equal we have 

This rule will be found convenient in practice to supply the 
defect in the number of arbitrary constants produced by equal 
roots. At present we are only concerned with their eflfect on the 
stability of the system. The terms which contain ^ as a factor 
will at first increase with <, but if m be negative, the term t^e^ 

can never be numerically greater than — . If w be very small 
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the initial increase of the terms may make the values of x and y 
become large, and the motion cannot be regarded as a small 
oscillation. But if the system be not so much disturbed that 



em 



is large, the tenns will ultimately disappear and the 

motion may be regarded as stable. If, however, the real parts of 
the equal roots are positive or zero, the terms will become large 
and the motion will be unatable. 

6. In some cases, however, the relations which exist between 
the coefficients are such that the terms which contain ^ as a factor 
are all zero. It is of some importance to discriminate these cases, 
for the stability of the system is then unaffected by the presence 
of equal roots. 

Let us suppose first that the determinantal equation has two 
roots only equal to m^, and let the terms depending on these be 

a? = ( JVj + N^) e'»l^ 

&c. = &c. 

Substituting in the equations of Art. (3) we have, following 
the same notation as before, 

j'jv;+5'jv;'+o'jv;"+... = o[ i. 

&c. = OJ 
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&c. = &c. 

To avoid entering more minutely than is necessary into the 
properties of linear equations, we shall assume that these equations 
for the given value of m lead to but one solution with two of the 
ITb arbitrary, unless the determinantal equation has more than 
two roots equal to w^. K in this unique solution the N^b are all 
zero we must have 

a'n,^bn; + .,.=^o\ III. 

&c. = 



6 EFFECT OF EQUAL ROOTS [CHAP. 

Since two of the constants N., N\ &c. are to be arbitrary, let 

N N' 
them be N., N' then since -^ = _ l ^e must have the minors a 

and 6 each equal to zero. Also since —^ = — ^, we shall have N" 

Cb C 

infinite unless c = 0. In the same way we may prove that all the 
other first minors are zero. 

And if the first minors are zero, we may show that two of the 
equations may be deduced from the others. Let the symbol 

A'jy\ represent the second minor, with the usual sign, formed 

)y omitting the rows and columns in which AB A'B occur. 

Then since the minors a, 6, c, &c. are zero, we have 



^U5'J-*-^''U'5''j+---=^ 



IV. 



&c. = 0. 

Omitting the first line of III. let us multiply the others by 
Afjy\ A"jy' ^^'t respectively. Adding the results, we have 

an identity. Hence the second equation may be deduced from 
the others which follow it. In the same way, the first equation 
may be deduced from the others. 

Rejecting the first two equations, let us transpose the arbitrary 
constants N. and N^' to the right-hand sides of the remaining 
equations. If there are to be only two arbitrary constants, these 
remaining equations must be independent; solying, we have 

with similar equations for the others. Hence the constaiits 
N^, N^', &c., are connected by equations of the form 

^.[IS]-^.'[i?]-^."[il]=«. 

SO that when any two are chosen as the arbitrary ones, the others 
may be deduced from them. 

If the determinantal equation has three roots equal to m^, and 
if the terms which contain ^ as a factor are all zero, the equa- 
tions III. must admit of a solution with three of the constants 
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JVj, N^\ &c. arbitrary. If these be i\^j, JV/, N^\ we see that the 
second minors 7y/r>i'/7/l4/iy must be zero. 

Next since 

* 
we mustliave 1 D/ry and I^/tv both zero or N"' infinite. 

Hence all the second minors are zero. And if the second 
minors are all zero, we have 



A" 



■ ABC 
A' BO' 
A"B'C' 



+ A 



Itt 



ABO 

a:fc 



+ ... = 0, 



and by similar reasoning three equations may be deduced from 
the remaining ones. We have then 



JV^. 



" BCD ' 


-^i 


' ACD ' 


+K 


' ABD ' 


-n;" 


' ABC ' 


BCD' 




A'C'U 




ABB 




A'ffC 


E'CB'. 




A" a'jy. 




A'S'iy, 




A"B"U'. 



=0, 



with similar equations. 

Since f[rn) is the determinant formed by eliminating N^, N^, 
&c. from III. we have 

df{m) ^df(m)dA ^ df(m)dB ^ ^^ ^ df(m)dA' ^ ^^ 
dm dA dm dB dm * dA' dm 



dm dm 



dm 



This vanishes when a, J, &c., a, &c. are all zero. If therefore the 
first minors of f(m) all vanish when ?» = m^, the equation /(m) = 

da 

dm 
first minors are zero. But 



has two roots equal to m^. In the same way ^ vanishes if all its 



d^/(m)_ d'A dadA 
dm dm dm dm 

vanishes if a, j— , &c. are all zero. If therefore the first and 

dm 

second minors of f{m) all vanish when m = m^, the equation 
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/(m) =0 has three roots equal to m^. It is evident the proposi- 
tion may be extended to any number of roots. 

The test that, when equal roots occur in the determinantal equa- 
tion, the terms in the valves of x, y , &c. which contain t as a factor 
sJiould he absent may he stated thus. If tJiere are two equal roots 
all the first minors must vanish. If three equal roots, all the first 
and second minors must vanish, and so on. In these cases the 
equal roots introduce merely a corresponding indeterminateness 

into Hie coefficients. 

• 

When there are more equal roots than there are rows in the 
determinantal equation, it is easy to see that there must be some 
terms in the integrals which contain ^ as a factor. 

[The following simple example will illustrate the application 
of this test. 

A particle is in equilibrium at the origin of co-ordinates under 
the action of forces whose force function IT is given by 

If the level surfaces are ellipsoids and the force acts inwards, 
it is clear that the equilibrium of the particle must always be 
stable. If then any equal roots occur in the determinantal 
equation, the test should show that the terms which contain t as 
a factor are absent. 

If T be the semi vis viva of the particle and if its mass be 
taken as unity, we have 

Omitting accents and forming the discriminant of —m^T+ U 
we have the following determinantal equation : 

A-m^ F E 

F B-m' D 



E D G-m 



s 



= 0. 



This is the "discriminating cubic" which determines the axes 
of the quadric Z7= c, where c is a constant. The conditions that 
two of its roots should be equal, i.e. that the quadric should be 
a spheroid, are well known to be 

. EF ^ FD ^ DE 

where m^ is equal to either root. These are just the conditions 
obtained by equating any first minor of the determinant to zero. 
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The conditions that three of the roots of the cubic should be 
equal, i.e. that the quadric should be a sphere, are 

These are the conditions that every second minor should vanish. 

In this example we have taken the case of a single particle. 
Similar remarks however apply when any system of bodies is 
disturbed from a state of stable equilibrium. The oscillations 
may be found by the method of Lagrange. The final determi- 
nantal equation may be conveniently formed by equating to zero 
the discriminant of ^m^T+ J7, where Tis the semi vis viva with 
the accents denoting differentiatioDS J^h regard to the time 
omitted, and U is the force function. ^tTs a known theorem that 
the existence of finite ^guaLroots does not afifect the stability of 
the equiliBrium77 Hence the conditions for equal roots must be 
such as to maSH all the minors equal to zero. Conversely, this 
theorem will often conveniently give the conditions that Lagrange's 
determinant has equal roots.] 

6. That there should be a diflference in the modes in which 
equal roots affect the motion is no more than we should expect 
a priori. Suppose the coefficients of the equation /(m) = to be 
functions of some quantity n, and that as n passes through the 
value Uq, two roots become equal to each other. Let the quadratic 
factor containing these roots be m* + 2am + /8, and let us consider 
only the case in which a and /8 are real. We have a* — ^=0 
when n = nQ. If a* — )8 change sign as n passes through the value 
n^ the roots will change from a trigonometrical to a purely ex- 
ponential form, which would indicate a change from oscillatory to 
non-oscillatory motion. The passage from one kind of motion to 
the other may be effected through a motion represented by expres- 
sions having the time as a factor. But if o? — fi does not change 
sign, for example, if it be a perfect square for all values of n, there 
will be no change from one kind of motion to the other, and in this 
case we should expect that the motion when the roots are equal 
will be represented by terms of the same character as before. 
Briefly, we may expect equal roots to introduce terms with ^ as a 
factor at the boundary between stability and instability ; and to 
introduce merely an indeterminateness into the coefiicients when 
the motion is stable on both sides. 

It is easy to show that in the first of these two cases the 
minors could not contain either of the factors of m*+2am + j8. 
For since a' — )8 changes sign, these factors are in one case ima- 
ginary ; and therefore if one factor occur in any minor the other 
must also be present. The minors would not only vanish, but 
must have equal roots also. But as in Art. (3), ^^f{^) = oil — ah 
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Hence if all the first minors have equal roots it is clear that either 
/(m) has more than two equal roots, or all the second minors must 
vanish. The latter is imJ)ossible unless f{m) has more than two 
equal roots. 

These general considerations are not meant to replace the 
proofs given in the last article, but merely to explain how a 
difference in the effects of the equal roots might arise. 

7. Summing up what precedes, we see that if a dynamical 
system have n co-ordinates its stability depends on the nature of 
the roots of a certain equation of the 27ith degree. 

If the roots of this equation are all unequal, the motion will 
be stable if the real roots and the real parts of the imaginary roots 
are all negative or zero, and unstable if any one is positive. If 
several roots are equal the motion will be stable if the real parts 
of those roots are negative and not very small, and unstable if 
they are negative and small, zero, or any positive quantity. But 
if, as often happens in dynamical problems, the terms which con- 
tain ^ as a factor Bxe absent from the solution, the condition of 
stability is tbat the real roots and the real parts of the imaginary 
roots of the subsidiary equation should be negative or zero. 

8. When the equation/ (2)) = is of low dimensions we may 
solve it or otherwise determine the nature of its roots; the 
stability or instability of the system will then become known. 
But if the degree of the equation be considerable this is not a very 
easy problem. We shall devote the two next Chapters to the 
consideration of two methods by either of which, without solving 
the equation, we can determine the conditions that the real roots 
and the real parts of the imaginary roots should be all negative. 
The determination of these conditions has, it appears, never before 
been accomplished.* The consideration of the equations of motion 
will then be resumed, and the form of the determinantal equation 
f{D) = when the forces admit of a force function will be^ more 
particularly investigated. 

* [These conditions for the cases of a biquadratic and a qnintic had been fonnd 
by the Author in 1873, and read before the London Mathematical Society in June, 
1871 See also the third edition of the Author's Rigid Dynamics, Art. 436.] 



CHAPTER II. 



Statement of the theory by which the necessary and sufficient tests of 
stability are found. Objections to this theory. Arts, 1 — 3. 

These tests shown to be integral functions of the coefficients. Art, 4. 

Method of finding these tests when the coefficients of the eqtiation are 
numerical, or when several terms a/re absent. Arts, 5, 6. 

All these tests shown to be derivable from one called the fundamental 
term. Arts, 7, 8. 

The fundamental termfotmd as an eliminant. Art, 9. 

A method offnding thsfundam^ental term by derivation from the funda- 
mffntaliemn4^ime degree lower. Art, 10. 

Another a/nd better method of doing the same by means of a differential 
equaiion. Arts, 11 — 13, 

1. The object of this Chapter has been explained at the end 
of Chapter I. Briefly, the criterion that the motion of a system 
of bodies should be stable is that the roots of a certain equation 
should have all their real parts negative. We propose to investi- 
gate these conditions. 

Let the equation to be considered be 

f{x) =Poi»" +pX"* + ,.. +Pn-i » + P»= 0. 
Let the real roots be a^, a,... and the imaginary roots be 

«! ± /8, V - 1, a, ± /3, V^=l, &c. 
Then 
/(^) =i>o (^- «i) (^ -«») ••• (^ - Sa^a; + ai' +^i')> &c. 

If then a^, a^, &c. a^, ofj, &c. are all negative, every term in 
each factor, and tnerefore in the product, must be positive. 

It is therefore necessary that every term in the equation 
f(x)=zO should have the same sign. It will be convenient to 
suppose this sign to be positive. 
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It is also clear on the same suppositions that none of the 
coefficients j9o,pj, ...^^ can be zero, except when the roots of the 

equation are all of the form + P^J— 1, or when some of the roots 
are zero. 

2. Let us now form the equation whose roots are the sums of 
the roots oif{x) taken two and two. Let this be 

n — 1 
where m = w — o~ • The real roots of this equation will be 

a^-k-a^^ a^ + a^, &c. 22^, 2afj, &c. and the imaginary roots will be 

c^i + flfj + ^1 V — 1, &c. It is clear from the same reasoning as 
before that if a^, a,, &c. a,, ofj, &c. are all negative, the coefficients 
P^ P,, &c. must all have the same sign. 

Conversely, if p^y p^... have all the same sign, the equation 

/{x) can have no real positive root, and if P^, Pj ... P^ have all 

the same sign the equation F(x) can have no positive root, and 

therefore f{x) can have no imaginary root with its real part 

positive. 

3. Ov/r first test of the stability of a dynamical system is thai 
all the coefficients of the dyna/mical equatiim f (D) = avd all the 
coejffhients of its derived equation F (D) = sh&uld have the same 
sign 

It should be noticed that though these conditions are all 
necessary and sufficient, they are not all independent. We obtain 
too many conditions. In many cases, however^ we can at once 
reduce them to the proper number of independent conditions, and 
when this is difficult we can have recourse to the second method, 
to be given in the next Chapter, which is free from this objection. 

In order to apply this method with success, it is necessary to 
have some convenient methods of calculating the coefficients 
P P P 

4. The first method which suggests itself is one similar to that 
usually given to determine the coefficients of the equation whose 
roots are the squares of the diflferences of the roots of any given 
equation. 

If 5j, iSj... be the sums of the first, second powers, &c. of the 
roots of the equation fx = 0, we have by Newton s theorem 

where p^ has been put equal to unity. If \, iS^... be the sums of 
the powers of the equation F{x) = 0, we have in the same way 



il] has positive coefficients. ]3 

we may also prove 

2, = (n -2)S, + 8', 

^,= {n-i)S,+ W,8„ 

^^ = (n-8)8, + 48,8, + S8*, 

2, = (n - 16) fi; + 58,8^ + 105,6; ; 

and the general relation can be found without difficulty. 

In this way we find 
P, = {n-l)p,. 

p_ («-l)(n-2)(n-3)(n-4) . (n-2)'(n-3) , 
* 1.2.3.4 P*'*' 1.2 P^P* 

But the process becomes longer and longer at every stage. 
We shall therefore proceed to point out some other methods of 
obtaining the coefficients. 

This method of proceeding has indeed been stated only because 
it proves in a convenient way that when p^ = 1, all the coefficients 
Pj, Pj, Pg... of the derived equation are integral rational func- 
tions of the coefficients jp^, p^... p^. 

5. The equation f (x) = being given, to calculate the coefficients 
o/F(x)=0. 

Put x=^y ±z and equate separately to zero the sums of the 
even and odd powers of z, we have 

/(y)+/"(y)^+r(y)g+-...-oj 

/'(3r)^+/"'(y)^+... =0 

Rejecting the root jg: = 0, let us eliminate «. Then the roots 
of the resulting equation in y are the arithmetic means of the 
roots of/ (a?) = 0. 
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If, on the other hand, we eliminate y we have an equation of 
an even degree to find z. This, putting 4a' = f, is the equation 
whose roots are the squares of the differences of the roots of the 
given equation. 

It may be thought that this elimination may prove tedious, 
but it will be presently shown that only the first and last terms of 
the result are really wanted. All the others may be omitted in 
the process of elimination, and thus the labour will be greatly 
lessened. The method is however most useful when the given 
equation has several of its terms absent. 

6. EaxLmple, To determine the condition that the roots of the 
biquadratic 

X* +px' + qx' + rx + s.= 
should indicate a stable motion. 

Applying the rule we have 

+ (2p"gr +pr + gr'- 4^) a?+ {p^ +p^r - 4ips) x+pqr-r^ - p's. 

The first four coefficients contain only positive terms, and need 
not be considered. If the last three coefficients be called P^, P^, 
Pg, we have 

p'P,- 4P,= (p^-2r)« + 2p*q + p\ 
pP,-iP, = {pq''2ry + p'r. 

If then Pg is positive, all the other coefficients are positive. 

The necessary and sufficient conditions of stability are there- 
fore that p, q, r, s should be finite and positive, arid 

P9=pqr''r^—p^s 
positive or zero. 

7. In forming the derived equation F(x) the only difficulty 
is to form the last term P^. For when this is knovm the other 
terms can he at once derived from it by an easy process. 

Let a, b, c... be the roots of /(a?) = with their signs changed, 
and let 

/(x) =py +p^x^' +....+ p^^. 
Let A stand for the operation 

da db dc 



.^ 
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Then since ~ = ahc ... we have obviously 
In the same way we have 



And generally 



A^2=i^ = iC^5=^ 



and so on up to 



Po i>, 







Po 



Let us now operate with A on any expression 

^(Po> Pi •••?»)> 

which has the same number of factors in every term. Let r be 
the number of factors, then ^ may be written 

"A ^ + (" ~ ^)^' ^, ■*" '" j *' 

Let P^ Pj ...P^ be the coefficients of the derived equation 
F{x), and let P^ = 1. Then since the roots of F{x) are the sums 
of the roots of /(a?) taken two and two, it is easy to see that 

&c. = &c. 

Thus when P„ is known, the other terms may be calculated 
without difficulty. The term P^ will be called the fundamental 
term of the equation. 

Example, Given in the case of a biquadratic 

•P«=PiP2P8-PoP.'-i'iX 
to calculate P^. 
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Performing the operation 

on Pj we find after division by 2, 

which is the result already given. 

8. It should be noticed that in the equation 

f{x) = ;)X +Pia;"'* + . . . +Pn = 0, 

if we regard a? as a number, p^p^ ... jp^ are all of equal dimensions. 
It follows from the theory of dimensions, that if any subject of 
operation be the sum of a number of terms of the form 

there must be the same number of factors in every term. For 
example, in every term of the expression for P„ we have 
OL + fi + y &c. the same. 

On the other hand, we may regard ^ as a quantity of one 
dimension, and in this case Po> J^i '^Pn ^^^® their dimensions in- 
dicated by their suffixes. We must therefore have ^+27+88+... 
as well as a 4-/8 + 7+ ... the same in every term. 

These two tests of the correctness of our processes will be 
found convenient. 

9. The whole derived equation being hnovm when the funda- 
mental term is knovm, it is required to find the fundamental term, 

[First Method.] 

If we write — a? for a; in any equation, we have a second equation 
whose roots are equal and opposite to those of the first equation. 
If we eliminate x between these two, we shall get a result which 
must be zero when the two equations have a common root. The 
eliminant must therefore contain as a factor the product of the 
sums of the roots of the given equation taken two and two. 

It will afterwards be shown that the last term P^ of the 
derived equation (whenp^ is put equal to unity) always contains 
the term 

with a coefficient which is positive and equal to unity. 
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Hence we have this rule, to find P„ . eliminate x between 

and divide the result by the coefficient o{pj)j)^...p^_^. 

It is obvious that the result may be written down as a 
determinant. On trial, however, it will be found more convenient 
to make the elimination by the method of eliminating the highest 
and lowest terms than to expand the determinant. 

10. Given the fundamental term of the equation derived from 

f{x)=p^'+p,^+... + p^_^0 (1). 

to find the fundamental term of the equation derived from 

p.x'^^p^x''-'^- .,.^p^,^x^p^ = (2). 

[Second Method.] 

Let Q„ be the product of the sums, two and two, of the roots 
of the equation (2) taken with their signs changed, so that Q^ is 
the same as the fundamental term of the derived equation and 
dififers from P^ only in having a suffix more convenient for our 
present purpose. 

Let Q^ be expanded in a series of powers of />^ : thus 

where <^<,, <^j, (f)^, &c., are all functions ofp^, p^, &c., which functions 
have to be found. Let the roots of f{x) = with their signs 
changed be a, J, c..., then 

C»-i=(a + 6)(a + c)... 

Let us introduce a new root, which, when its sign is changed, 
we shall call r. Then 

C„=(a + Z>)(a + c) ... (r+a)(r + 5)... 
= %■' (Par""" +P^r'"^ + . . . + 2>n-i). 

Po 

This value of Q^ must be the same as that given by the series 
when we write 

Po* Pi + ^oy P2+'^Pi> &c., ;>,_!+ rp._,, rp,_^, 

respectively, for p^, p^, p^, &c., p„_^, p^. Equating the coefficients 
of the terms independent of r we have 

R. A. 2 
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Equating the terms containing the first power of r we have 

Substituting for <f>^ we have, by a known theorem in the Diflferential 
Calcuhis, 

Equating the terms containing the second power of r we have 

d d 



Fa 



and so on. Thus we have 

^^ ^""' Po ^""V' dp^ ^'dp^ ^^dp^J p. 

If we examine this process, we see that when Q„_j is known, 
we may at once write down the terms independent of p^ and the 
coefficient of p^. The process to find the coefficient of p^ is 
longer, but it may be much shortened by the Consideration that 
when Pq=1 the result must be an integral function of the coeffi- 
cients. We may therefore omit all terms as soon as they make 
their appearance, which do not contain the factor jp*„_i, for we 
know that such terms must disappear from the result. 

This method is not so convenient as the one which will be 
presently given to find the coefficients of the higher powers of p^. 
But it is useful as showing that Q^ contains the term 

PJ>i"'Pn-l 

with a positive integral coefficient equal to unity. This will be 
clear from the consideration that the term independent of p„ in 

Q^ is obtained from Q^.^ by multiplying by -^^^ . If therefore Q^ , 

contains ^^^S^^', Q. must contain the term P^P^"'P--^ . No 

Po . . Po 

other term can be formed which is equal to this with an opposite 

rign, for the terms which enter by the other processes to be per- 
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formed on Q^_^ all contain jp» as a factor. Now Q, = - ; therefore 
Q^ contains the term ^-^' , Q contains ^^2?/-?, and so on. 

Po Po 

It has been shown that every term in p^^Q^ has the same 
number of factors (Art. 8). It follows from this reasoning that 
this number of factors is n — 1. 

11. To find the fundamental term of the derived equation by 
means of a differential equation. 

[Third Method.] 
The fundamental term required is a factor of the eliminant of 

Let »* = y, then we have 



J-i ?-t 



even, 



Pif '+py"'+ =o[ 

w— 1 w— 8 > 

• i'.y'+i'.y'+ =o' 

If we write y + dy for y the result of the elimination must be 
the same. Hence if we make 

dy 2^0' dy \2 / ^*' L 

dy [2 ^JP^' dy-[2 ^)P"'^\ 

dy 2 P"' dy~ 2 P" ' 
dp^_ n-l t7j?,_ n-3 _ ' 

'dy~ 2 P" dy 2~^»' 

and if ^ be the eliminant, we have 

dy 

It follows that whether n be even or odd, ^must satisfy the 
equation 

2—2 



n odd, 
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dE dE 



-/ dE dE\ 

^-/ dE ^ dE\ 



+ 



^ = 0. 



We may make the elimination by multiplying the two equa- 
tions by y, y*..., until we have as many equations as we have 
powers to eliminate. If in the determinant thus formed, we 
multiply out the terms in the diagonal joining the right-hand top 



■S-i ■?+! 



corner to the left bottom comer, we get when n is even p * p * 

n-l n-1 •* * •* * 

and when n is odd p^ ' i'o " • -^^^^ Qn ^^st contain w — 1 factors 
and be of the n — ^ — th degree. Hence when n is even E=- cp^Q^ 
and when n is odd E= c^„, where c is some constant. 

Now -T— does not occur in the above dififerential equation. 

Hence treatmg p^ as a constant, we see that Q^ must satisfy the 
differential equation 

V^- dp, + ^-» d'pj + (,^"- c?p^ + ^- rfp^ J 



^3(,..^^,^.^^ + &c. = 0. 



12. We may show that p^'^Q^ is a symmetrical function of the 
coefficients p„ p,..., p, ancfthe same coefficients read backwards. 
Let a, 6, c... be the roots of f(x)=0 with their signs changed, 
then <?n = (« + ^) (a + c).... 

If now we read the coefficients in the opposite order, the roots 
of the equation thus formed will, when their signs are changed, be 

- , ^.... If Qfn be the fundamental term of the equation derived 

from this, we have 



'^H^Dil^D-- 



Since — = abc... we see that 



i 
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We may therefore infer thsitp^'^Q^ also satisfies the differential 
equation 

/ dE ^ dE\ ^ ^f dE ^ ^^\ . o / ^^ . d^\ . jp A 

[p^dJ^^^^)Hp^w^'''^^^^ 

13. We may use either of these differential equations to find 
Q„ when Q^^ is given. 

Let the first differential equation be represented by 

and let 

Q^ = A^ + Aj)^ + Aj>^^ + ... , 

where A^^ A^, ... are functions ofp^, p^^ &c. The value of A^ has 
been proved in Art. 10 to be 



A = <2«-. 



V 



n-l 



Pi 



To find the other coefficients of the powers of ^„ substitute 
this value of Q„ in the differential equation ; we have 

= VA, +PnVA +i>»-3^i 

+ &c. 
Equating the several powers of ^^ to zero, we find 



A= 


1 


t 


2A,= 


1 

P»-i 


1 


3A = 


1 


vA> 



&c. = &c. 

Thus by one regular and easy process each term may be 
derived from the other. 

In performing this process we may omit every term in the 
subject of operation which does not contain p^.^. For p^^ can be 

introduced only by performing -r— , and since p» is absent from the 

coefiicients, this operation yields nothing. 
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Ill this way we find 

P'Q»'=PiP»-p»P,' 

P*Q, = PtPiP»-PoP,-PiP*> 

P»Qi '=PiPtPJ>*-P»PtP4 -PiP* 

-P, {-P„P^» +PiP* - ^PoPxPd 



+ f^i- 2p/). 



To illustrate this process, consider how Q^ is obtained from Q^. 
The first line is formed by multiplying the line above by p^, this 
is A^ To find the coefficient of —jp^ we operate with 

{ps^+p4)'-Kp^i-^P'^) 

on such of the terms in the line above as contain p^ and then 
divide by p^. Performing the same operation on the coefficient of 

("JPt) we obviously obtain the coefficient of r^^ . 

In M. Serret's Cours cPAlghbre Sup^r{eure,Noie III., there will 
be found a method of forming the last term of the equation to the 
squares of the diflFerences, which suggested the method used in 
Art. 13, of substituting in a differential equation, if only a differ- 
ential equation could be found. [See also Dr Salmon's Higher 
Algebra, Arts. 60, 64, and 72.] 



CHAPTER III. 



Statement of the tlieory by which the conditions of stc^Uity of a dynamical 
system t(nth n co-ordinates a/re made to depend on 2a condition's. 
Arts. 1—4. 

A rule by which these conditions may be derived one from, oimther^ to- 
gether with certain other true but not independent conditions. Arts, 

A rule by which, when the coefficients of the dynamical equation a/re letters, 
the 2q conditions of stability may be inferred, one from another, 
without writing down any other conditions. Arts. 9, 10, 

A meiliod by which certain extraneous factors m,ay be discovered and 
omitted. Arts, 11, 12. 

Consideration of the reserved case in which the dynamical equation has 
equaJi and opposite roots. Arts, 13 — 18. 

A few geometrical illustrations not necessary to the argument. Arts, 
19—25. 

Application to a Dynamical Problem, Art. 26. 

1. It has been shown in the first Chapter that the stability of 
a dynamical system with n co-ordinates oscillating about a state of 
steady motion depends on the nature of the roots of a certain 
equation of the 2n^ degree which we may call 

The system is stable if the real roots and the real parts of the 
imaginary roots are all negative. Now Cauchy has given the fol- 
lowing theorem of which we shall make some use. 



Let z^x-^-ysl — 1 be any root, and let us regard x and y as 
co-ordinates of a point referred to rectangular axes. Substitute 
for z and let 

f{z)^P+Q^^^. 

Let any point whose co-ordinates are such that P and Q both 
vanish be called a radical point. Describe any contour, and let 
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a point move round this contour in the positive direction and 

p 
notice how often j^ passes through the value zero and changes its 

sign. Suppose it changes a times from -f to — and /8 times from 
— to +. Then Cauchy asserts that the number of radical points 
within the contour is ^ (a — )8). It is however necessary that no 
radical point should lie on the contour. 

2. Let us choose as our contour the infinite semicircle which 
bounds space on the positive side of the 4xis of y. Let us first 
travel from y = — oo toy = + oo along the circumference. 






we have changing to polar co-ordinates 

/(2?)=j>/*(cosw5 + 8in?i^V— 1) -I-... 
Hence 

Q = py sin n^ + p,r"-* sin (n- 1) ^ + ... ) " 
In the limit, since r is infinite, 

g = cot nd ; 

-7:z vanishes when n6 = (2/c + 1) — , i,e. 

^=±-?. ±-%, ±-l (A); 

pr is infinite when n^ = 2/c ^ , i.e. 

^-^' ^n'l^ ±n2' ±n2 ^^^' 

The values of in series (B) it will be noticed separate those 
in series (A). 

P . 

When is small and very little greater than zero, jr is posi- 

tive, and therefore changes sign from + to — at every one of the 
values of in series (A). If n be even there will be n changes 
of sign. If n be odd there will be w — 1 changes excluding 

0= ±-^ y in this case ^ is positive when ^ is a little less than -^ , 
and negative when 5 is a little greater than -^, 
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Let us now travel along the axis of y still in the positive di- 
rection, viz. from y^ + oo toy= — oo. Since ^ = it will be 
more convenient to use Cartesian co-ordinates, we have, since 



and 



and z^ysf —I, 

The condition that there should be no radical point within the 
contour is that this expression should change sign through zero 
from — to + as often as it before changed sign from + to — on 
travelling round the semicircle. If n be even the numerator has 
one more term than the denominator, and when 'p^ and 'p^ have 

P . 

the same sign, -^ begins when y is very great by being negative. 

In order that it should change sign through zero n times, it is 
necessary and sufficient that both the equations 

Pn -Pn-^f -^Pn-iSt " • • • = 0, 
Pn^xV - Pn-^y"" +jPn-62/' " • • • = 0, 

should have their roots real, and that the roots of the latter should 
separate the roots of the former. 

If n be odd, the numerator and denominator have the same 

p 
number of terms, and when p^ and p^ have the same sign, jz begins 

when y is very great by being positive. In order that it should 
change sign through zero from — to -h n — 1 times, it is necessary 
and sufficient that the same two equations as before should 
have their roots real, and that the roots of the former should 
separate the roots of the latter. 

In order then to express the necessary and sufficient conditions, 
that f(^z) = iruiy hive no radical point on the positive side of 

the axis of y, put z = y v — 1 and equate to zero separately the 
real and imaginary paHs. Of the two equations thus formed, the 
roots of the one of lower dimensions must separate the roots of the 
other. It is also necessary that the coefficients of the two highest 
powers of z in f (z) should have the same sign. 

3. It has been stated that p^ and p^ the coefficieilts of the two 
highest powers in f{z) must have the same sign. It is easy to see 
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P 

that, if they had opposite signs, ^ would change sign through zero 

2» times as we travel round the contour. All the radical points 
of the equation would then lie on the positive instead of the 
negative side of the axis of y. 

It has also been assumed that no radical point lies on the 
contour. It has therefore been assumed that/(«) =0 has no root 
of the form z = y V— 1. It will be more convenient to consider 
this exception a little further on. 

4. It is required to express in an analytical form the conditions 
that the roots of an equation f, (x) = rrvay he all real, and may 
separate the roots of another equation f^ (x) = q/* one degree higher 
dimensions. 

To effect this, let us use Sturm's theorem reversed. Perform 
the process of finding the greatest common measure of f^ (x) and 
f (x), changing the sign of each remainder as it is obtained. Let 
the series of modified remainders thus obtained he f^(x), f{x)^ 
' &c. Then it may be shown that when any one of these functions 
vanishes, the two on each side have opposite signs. It is also clear 
that no two successive functions can vanish unless f (x) and f^ (x) 
have a common factor. This exception will be considered pre- 
sently. 

Hence in passing from a: = — oo to + oo no variation of sign 
can be lost except when f (x) vanishes. If a yatiation is lost it is 
regained when x has the next greatest value which makes f (x)^ 
vanish unless ^ (a?) = has a root between these two successive' 
roots of ^ (x) = 0. Hence this rule : — 

The roots of the equations fj(x) = 0, fj(x) = 0, vrill be all real 
and the roots of the latter will separate those of the former, if in 
the series 

^i(x), f2(x), faCx)... 

as many variations of sign are lost in passing from x = — oo to 
X = + 00 cw there are umts in the degree of the equation f^ (x) = 0. 

We have supposed the variations of sign to be lost instead of 
gained in passing from a; = — oo to + oo . That this may be the 
case the signs of the highest powers of f (x) and f^ (x) must be 
the same. 

These functions are alternately of an even and odd degree, the 
condition that the whole number of variations of sign may be lost 
in passing from a; = — oo to +oo may be more conveniently ex- 
pressed thus : — The coefficients of the highest powers of x in the 
series 

f,(x), f,(x), {,{x) ... 

must all have the same sign. 
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5. The process of finding the greatest common measure of 
two algebraic expressions is usually rather Ions:. We may in our 
ca^e shorten it materially by omitting the quotients and perform- 
ing the division in the following manner. Let 

Sx (^) =i>oaj" -V^'^ + A^"^ - • • • 

/, (^) ^v.^""" -i>8^""' +i>5^"^ - • • • 

then, since p^ is positive, it easily follows by division that 
where 'A^p^'p^-'p^'p^, 

&c. = &c., 

so that by remembering this simple cross-multiplication we may 
write down the value of f, (x) vritlwut any other process than what 
may be performed by simple inspection. In the same way ^(a?), 
&c. may all be written down. 

6. Ex. 1. Express the conditions that the real roots and 
real paxts of the imaginaiy roots of the cubic 

ic' -hpx* + ja; 4- r = 
may be all negative. 

f {x) = a?' - qx, 

f^{x)^p3?-r, 

fs{^)='{Pi-r)x, 
f,(x)={pq-r)r. 

The necessary conditions are that 

p, pq—r and r 
must all be positive. 

Ex. 2. Express the corresponding conditions for the bi- 
quadratic 

aj* + pa;' -hqx^ + rx + s = 0, 

fi{x)=x^ —qa^ + s, 

f,(x) = {pq-r)a^ -ps, 

A{^)={{pq-r)r'-p's}x, 

A{^) = {(pq'-'r)r^p^ps. 
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The conditions are that 

p, pq — r, {pq — r) r—p^a and s 
must be all positive. 

These are evidently equivalent to the five conditions that 

P, J, r, 8, {pq-r)r-p\ 
should be all positive. 

In both these examples all the numerical work has been 
exhibited. 

7. Since the coeflBcients of the highest powers of x infjjv) and 
/j {x) are p^ and p. we see that the condition that p^ and p^ should 
have the same sign is included in the general statement that all 
the coefficients of the highest powers should have the same sign. 
If the function /(a;) be of n dimensions we thus obtain n necessary 
and sufficient conditions. 

On examining these conditions in the cases of the cubic and 
biquadratic it will be seen that they cannot be satisfied if any one 
of the coefficients of the given equation should be negative. 

8. Although the theorem in its present form gives n conditions 
as the proper number for an equation of the n^ degree, yet it is 
important to notice that it gives other conditions also which are 
true and may be useful. It has been shown in the second Chapter 
that all the coefficients of the equation f(x) = must be positive, 
hence the roots oif^{x) = must all be positive. It may be shown 
also, that the roots of each of the functions fAx), f^{x), &c. are 
separated by the roots of the function next below it in order. 
Hence the roots of all these functions must be positive, and there- 
fore in every one of the functions the coefficients of all the powers 
must be alternately positive and negative and not one can 
vanish. If however f^{x) and f^{x) have one or more common 
factors some of the functions /3(a?),/^(ir), &c. will VhoUy vanish. 

9. When the degree of the equation is very considerable 
there is some labour in the application of the rule given in Art. 5. 
The objection is that we only want the terms in the first column, 
and to obtain these we have to write down all the other columns. 
We shall now investigate a method of obtaining each term in the 

first column from the one above it without the necessity of tvriting 
down any expression except the one required. 

We notice that each function is obtained from the one above it 
by the same process. Now 

fsi^) = (P1P2 -P0P3) ^"^ - (PiP* -i^oPfi)^""* + • • • 



Po. 


Pi> 


the values 




Pv 


PiPi -PoP<^ 
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The first and second lines will be changed into the second and 
third lines by writing for 

P2> P^ &^- 

i^3> PiP^-P^6>&^' 

If then in any term of any function we make these changes, we 
obtain the corresponding term of the function next in order. 

10. Example. Express the conditions of stability for the 
quintic 

f{x) = p,x^ + p^x* + p^x^ + jPjX' + p^x + p,. 

We have 

♦ 

/i(^)=i>o^ + --'- 

f,{os)=p^x*+... 

f3(^) = (PlP2-PoP^)^'^ 

A(^) = {(PiP^-PoPM-PiiPiPi-PoPs)]^' 

M^)=[{(j>J>2-PoPs)P3-Pl(PlP4-PoP^}(PiPi''PoP^ 

-{pxP^-PoP^y^p^^y 

fM = \.{{PiP^ - PoP^P^ -P. iP^P^-PoP^] {PiP, -PoPd 

11. On examining the conditions as given in the cases of a 
biquadratic and quintic it will be apparent that several contain 
the previous conditions as factors. Thus the analytical expressions 
are rendered much longer than is necessary. It is now proposed 
to investigate a method of discovering and omitting these extraneous 

factors as they occur, and thus obtaining the required conditions in 
their simplest forms. 

Let the coefficients of the several powers of x in the functions 
be when taken positively 



fl{0')=P<» 


Pv 


p,. 


Pe •*' 


M'«)=Pv 


Pz" 


p» 


Pi'*' 


/.(<«) =^, 


A', 


A", 


,/L • . . 


M^)=B, 


b; 


b; 


jy// 



&c. = &c. 
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Let us first fiod which of these terms contain pi as a factor. 
Patting />i = and using the rule in Art. 5, the series become 



Po. 


?.. 


Pi' 


p., ... . 


0, 


p.. 


p»> 


p, 


PoPt> 


— PoPsi 


-PoPk 


-PoP» ••• • 


Popi. 


-PoPzPt' 


-PoPsP7> 


-Pi>pipt> — • 


0. 


0, 


0, 


0, .•.. . 


0, 


0, 


0, 

&c. 


0, ... . 



Hence the C's and D's ail vanish and therefore contain p^ as 
a factor. By the rule in Art. 5, the E^s contain p^, the F's con- 
tain p^^, the GPspi, the ITsp^, and so on. 

But sinc6 each line is formed from the preceding by a uniform 
rule, it follows that the D*a and E's contain A as a. factor, the 
F'a contain A\ the O's contain A^, the ITs contain A\ and so on. 

The factor A in the i)'s and E's has its origin in the factor p^ 
which occurs in the (7s and D's and would not appear if that 
factor had been omitted when the C's and D's were formed. The 
factor A in the D's and Es in the same way gives rise to the 
factor B in the -E"s and Es. So that if we take care each time 
we perform the process described in Art. 9 to omit the common 
factor p^ whenever it occurs, all these subsequent factors will never 
make their appearance. 

We shall now show that if these factors are omitted, the 

dimensions of the n*** function j^ (a?) will be n — ^ — . First consider 

the actual dimensions of each function before the factors are 
omitted. If we examine the rule by which each function is 
derived from the preceding, it will become evident that, the 
dimensions of each letter being indicated by its sufiix, the 
dimensions of any function are equal to the sum of the two pre* 
ceding + 2. 

In the following table the first column indicates the function. 

In the second column will be found the dimensions of the leading 

coefficient of that function when calculated by the rule in Art. 5. 

In the third column will be found the dimensions as given by the 

w— 1 
formula n — ^ . In the remaining columns are the dimensions 

of the extraneous factors jf?^, -4, B, &c. introduced into each term. 



\ 
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fM 
















M^) 


1 


1 










fM 


3 


3 










fM 


6 


6 










M^) 


11 


10 


1 








fM 


19 


15 


1 


3 






fM 


32 


21 


2 


3 6 






fA^) 


53 


28 


3 


6 6 


10 




fM 


87 


36 


5 


9 12 


10 


15 


&c. 


&c. 


&c. 


&c. 


&c. &c. 


&c. 


&c. &c. 



Each term in the second column is the sum of the two terms 
just above it +2. The n**" term in the third column is equal to 
the term just above it + (n - 1). In all the other columns each 
term is the sum of the two terms just above it. The last term in 
the n*** row is equal to the {n—Z)^ term in the third column. 
We wish to show that any term in the second column is equal to 
the sum of all the terms in the same row to the right of that term. 
It is not diflScult to show from the data just given that if this be 
true for any two adjacent rows, it is true for all the others, and 
hence we may assume it to be always true. 

It is clear that these extraneous factors may be omitted since 
by the conditions already expressed they are all positive. When 
omitted as they occur, the dimensions of the rC^ function has just 

been shown to be n — — . It is easy to see that the conditions 

thus reduced must contain the terms 

Pv PiP2y PiP2Psy PiPiPzPi^^^' 
Now if we take any one of these as 

PiP2PsPiPf^Pe (1), 

and operate by the rule in Art. 9, we have 

iPiP^ - PoP^Ps (PiPi - PoP5)P, (PJ>6 - PoP,)Pv 
which contains the term 



PiP2PsPiP4P^PiPePi (2). 

Thus we have jp, introduced as often as there is a factor in (1) 
with an odd sufBx. But it should be introduced only once. These 
extra p^'s are the extraneous factors to be omitted. Each of these, 
if left, would appear as the factor jOj^^ ^PoPs i^ ^^^ ^^^^ condition, 
and be still more complicated in the next after that. 



/ 
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In order then to obtain the several conditions in their simplest 
form it is only necessary after performing the operation described in 
Art 5 or Art. 9 to divide by /?,*, where k is one less than the number 
of fajctors with odd suffixes in the condition operated on. 

12. Example. Express the conditions of stability for the 
sextic 

We have 
f^ («) = iVJPJP^ - V^Vz - VxV^ +i>o;>iP5) a?' + . . . 

/ft (a?) = {VxVJPzVi. - Ai^8>4 -'PxVl + ^PoPrPiP, -PiP^P, +l>oP,PzP5 

-PoP^ + PiP^P^ - PoPdP^P^] «* + ... 

/e (^) = [P1P2P8P4P5 -P^P^P^P^ - PiPaP^ + "^PoPrPiP^ - P1P2 V 

+ PoP,PJ>^* - i>oV + ^PiPJ>5Pe - 3jPoi>i;>8;>6P6 -PiPiPsPe 

-^PoPsPe+PiPsPiP, -PiPe] ^ + &C. 

/y (a?) = coefficient of a? in /g a? x by p^. 

13. In the preceding theory two reservations have been made. 

1. In applying Cauchy's theoreni it has been assumed that 
there were no radical points on the axis of y. 

2. It has been assumed that P and Q have no common 
factor, so that none of the functions /j,/,, &c. vanish absolutely. 

If any radical point lie on the axis of y, it is clear that f{z) = 
must have a.factor of the form (z^ + a')*". Let f{z) = (z^ + a^*" ^(«). 

In this case when we put « = y V — 1, we have 

/(«) = (a'-3^y(P'+<2'V^); 

Thus P and Q have a common factor, and we are warned of the 
possible existence of radical points on the contour by the total 
vanishing of some one of the fu notions /^jj^j/j, &c. 

The two reserved cases may therefore be included in one. If 
f{z) = be the equation furnished by dynamical considerations, 
we form the functions /j, f^, f, &c. If all these be finite, the 
question of the stability of the system has been answered. If any 
one vanish absolutely, f (x) and f (x) have a common measure, 
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and we must add some further considerations. It will be con- 
venient to examine separately the dynamical efifect of the roots 
which do not and which do enter through the greatest common 
measure. Let us begin with the former. 

14. Following the same notation as before, we have 

P^±A{y)=Pn-Pn^+ 1 



If then/(2j) have two roots, viz. +(A + & V — 1), which are equal 
and opposite, then f^ (y) and /, {y) must have two common roots, 

viz. + — . The common measure therefore of /, (y) and 

~ V — 1 

f^ (jf) contains all the roots of f(y V — 1) which are equal and 
opposite. Conversely the greatest common measure of P and Q is 
necessarily an even function of y, and if it be equated to zero, its 
rools are necessarily equal and opposite. These roots must also 

satisfy /(yV^) = 0. 

Let this greatest common measure be -^ (y") = 0, and let y*" be 
the highest power which enters into it. Also let 

then ^ (-s^) is a function which, as has just been shown, has not got 
two roots equal and opposite, and to this function we may apply 

Cauchy's theorem without fear of failure. Putting « = y V — 1, let 

P' 

Then we wish to express the condition that jy should change 

sign from — to + through Jero t* - 2r times if n be even and 
w — 2r — 1 times if n be odd. But 



/(2/V-i)=P+(2V^, 

and /(y V^^i) =^{y'){P' + Q' V":^); 

P 

Thus the number of changes of sign in -^ is exactly the same 

p 

as that of -q . The factor -^ (y*) will run through all the functions 

/a (y)> ft iy)> &c- obtained from /^ {y) by a process which is equiva- 
lent to that of finding the greatest common measure of f^ (y) and 

R.A. 3 
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/ (y). The changes of sign of this factor will therefore not affect 
the number of variations of sign in the series /i> /,, yi* ^^ 

The last factor which is not zero is /^j.,, (y) if « be the dimen- 
sions of/i (y). 

^ence if we omit the considerations of the vanishing factors and 
apply the same rule cw before to the n + l — ^r remaining fa/ytors^ we 
can express the condition that the proper nuTnher of changes of sign 
from — f + have been lost through zero in the function^ (z),t.e. that 
the roots not given by the vanishing of i^^^^ are all of the character 
to ensure stability, 

15. Let us next consider the effect on stability of the roots 
indicated by the absolute vanishing of one of the subsidiary func- 
tions. This function must be of the form 



where n is even. The corresponding factor of f{z) is 

^ W = ?n + ?tMI«' + ?»^«* + — • 

The roots of this equation are two and two equal with opposite 
signs, it is therefore necessary for stability that no root should 
have any real part. To express this condition, draw a straight 
line parallel to the axis of y at an indefinitely short distance from 
it, viz. x^h. Let us apply, in the same manner as before, Cauchy's 
theorem to the contour formed by this straight line and the infinite 
semicircle on its positive side. Putting <^ = A + is?', we have 

F {z) = y. + 2?„>/ + y^/^ + ^qnj^z'^ + . . . 

The two functions are therefore, omitting the positive factor A, 

e = 2j,.^-4j„.y+&c. 

Now/- (y) and/j (y) are what P and Q become^ when arranged 
in descending powers of y and the coefficients of their highest 
powers Tnade to have the same sign. Hence 

The rule described in Art. 4 will now become the same as 
that usually called Sturm's theorem. We are to seek the greatest 
common measure of/ (t/) and its differential coefficient, and make 
the coefficients of the highest powers of y in these two and in the 
series of modified remainders all positive. 
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That we should have been led to Sturm*s theorem in this ease 
is just what we might hare expected. For to express the condi- 
tions that the roots of 

are all of the form ± r V — 1 is the same thing as to express the 
conditions that the roots of 

are all real. 

16. There is however another mode of proceeding. Suppose 
we have calculated the functions /j, /j, &c. for the general equa- 
tion 

and find when the values of p^, p^, &c. are substituted that some 
one function say j^ of the series absolutely vanishes, and therefore 
also all the functions which follow it. Then operate on each of 
these vanishing functions with 

repeating the operation until we obtain a result which is not zero. 
If we now replace these vanishing functions by these results we 
may apply the rule of Art. 4, just as if these were the functions 
supplied by the process of the greatest common measure. As this 
process is not so convenient as that already given it is unnecessary 
to consider it in detail*. 

17. As a numerical example, let us examine whether the roots 
of 

/ (x) = ar* + 2x' + 4aj' + 4a?' + 6x* + 6aj' + 7a?' + 4a; + 2 = 

satisfy the conditions of stability. In order to show the working 
of the method it will be necessary to exhibit all the numerical 
calculations. We have by Art. 5, 

* The function /^ {x) vanishes because the equation /(2) = has two roots equal 
and opposite. If we put z=s^ + hy where h is as smaU as we please, this ^peculiajity 

will disappear. Thus if the values of z are of the form ^{a^fi J -^) ^^ corre- 
sponding values of z' are -h^a^p »J -1, These values of e will indicate stability 
if a be zero and instability if a have any value positive or negative. If ^ be as smaU 
as we please and positive, the values of z' will indicate stability or instability under 
the same circumstances. We may therefore apply the rule of Art. 4 to the func- 
tion f(^+h) instead of /{z), provided we retain only the lowest powers of h which 
occur. Hence all the functions fi(x), fq{x)...fr^i{x) which do not vanish are 
unaltered. To find what function will replace fr {x) we must increase by /i all the 
roots of /(2)=s0 when their signs have been ehuiged. This may be effected by 
performing on / (z) the operation represented by A in Art. 7. The rule in the text 
therefore foUows from the one given in that article. 

3—2 



\ 
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/, (x) = 4a?* -ex' + lOa?- 4, 
/ (a?) = 4a:^ - 4a;' + 8a?, 
/.(a?) = 8aj*-8aj' + 16. 

Here we find/e (a?) to be absolutely zero, accordingly by Art. 15 
we replace it by the differential coefficient of f^{x), this being 
Sturm's rule. We have therefore 

/,(a?)=8(4a^-24 

/, (a?) = - 8M20. 

We see that the two last of the coefficients of the highest 
powers are negative. The roots therefore do not satisfy the con- 
ditions of stability. 



As another example, take the equation 

j^U afl + ai'+6x^ + 5a?-^ 11a? + 6a? + 6. 

Here 

X(a?) = a;«-6aj*+lla»-6, 

/, (a?)=aj^-5a?"+6a?, 

/,(x) = a;*-6a?+6, 

/,(a?) = 0. 

Replacing/^ (a?) by the differential coefficient o{f^(x), we have 

/4(aj) = 4a?»-10a?, 

/, (a?) = 10a;' - 24a?, ! 

/^ (a?) = 4a?, 1 

/,(a?)=4.24. 

Here all the coefficients of the highest powers are positive, 
hence the roots satisfy the conditions of stability. 

It is clear that when the coefficients are numerical the rule 
given in Art. 5 is the most convenient, but when the coefficients 
are letters, the rule in Art. 9 will be found preferable. 

The process would be simplified by omitting the alternate posi- 
tive and negative signs of the terms in each line. 
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18. It may be interesting to express the two subsidiary funo- 
tions fj (x) and f, (x) in terms of the roots of the given equation. 

Let ttj, ttj, a,. . .a^ be the roots of the given equation/(a?) = 0, a[o 
that 

f(x) = (a?— a,) (a; — aj(a; — a,)... 

Then it is evident that 
±f {x "J^l) == aj* -^pjxT^ +P4p^'^ + . . • 

= J {x + a^)(x + a^) ... + J(aj-aJ {x-^a^ ... 

It may be shown* that 

_ ^ K + aJK + a,)...(a, + aJ . . . v 

19. The foUowing propositions are not necessary to the main 
armiment, but as they illustrate geometrically the propositions in 
this chapter it has been considered proper to state them very 
briefly. The demonstrations will therefore be much curtailed. 

The equation being 

we put as in Art. 1, *= a; + y V— 1, and thus obtain two curves, 
whose equations expressed in polar co-ordinates are 

P =i?o^ cos n^ +;)jr*~' cos (» - 1) ^ + . . . = 0| 
Q = jPor* sin w^ -\-'p;f'''' sin (n - 1) ^ + . . . = Oj * 

These intersect in the radical points of the equation/ {£) = 0. 

20. If we trace these curves we find that the curve P = 
has n asymptotes whose directions are given by cos nO = 0, ie. 

p_ 1 TT 3 TT 5 TT 

' w2' n2' n2' 



* Let ns assume 

Pia^^+p^it^^ + &Q.=Ai {x- a^,..{x- a^ + A^ (oc- ai)...(«- oj + &o., 

where J^, A^, &q, are constants whose values have to be found. Putting opso^, 
we have JPi^i*"^ +1^8^*'' +&o.=Ai{ai-a^... {oi - a J. 

But since as*+2)ia;*"i + ... = (a6-aj (a8-a8)...(«-aj, 

we have by putting x=ai and a;= - o^ 

«i*+l'i^'*"^ + — = 0, 

Subtracting the second of these results from the first we find J^ to have the value 
given in the text. 
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These asymptotes all pass through the same point on the axis 
of X, viz. 05 = — -2l , It is also clear that only one branch of the 

curve can go to each end of an asymptote. Similar remarks apply 
to the curve ^ = 0, the directions of its asymptotes being given by 
sinn^ = 0, i.e. 

*^" ' w2' n2' w2' 

21. From these simple propositions we might, if it were 
worth while, deduce that every equation must have a root. The 
asymptotes of Ihe two curves P = 0, Q = are alternate, and no 
two branches of the same curve can approach the same end of an 
asymptote. By sketching a figure, it may be easily shown that 
some branch of the P curve must cut some branch of the Q curve. 

22,^ Let us next consider the intersections of the curves P = 0, 
Q = 0. 

If we transform the origin to h, k, we put x = h + ^, y = k + i]. 

This is the same as expanding f(h + iV— l + f + iyV— 1), and 
collecting into two parcels the real and imaginary terms. Let the 
expansion be 

A+A(f+^^^^) + A(f+'?V-i)>+... 

where A^A^... are of the form 

c (cos a + sin a V^^). 

If we put f + ^7 V— 1 = r (cos ^ + sin ^ V— 1), we have 

P = CoCos a^ + c^r cos {0 + a^) + cy cos (25 + a,) + ...) 
C =" Co sin a^, + c^r sin (5 + aj + c^r^ sin (2 5 + a^ + ... J * 

If the point {h, h) be a point of intersection c^ = 0. If the inter- 
section be a double point on either curve, the terms of the first 
degree must be zero, therefore c^ = 0, and the origin is therefore 
a double point on the other curve also. 

It is not difficult to show that if the intersection be a multiple 
point of any degree of multiplicity on one curve, it is a point of 
the same degree of multiplicity on the other curve. The tangents 
to these branches all make equal angles with each other, the 
tangents to the P and Q curves being alternate as we travel round 
the point of intersection, If the intersection be not a multiple 
point on either curve, the branches cut at right angles. 

Let us travel round a point jof intersection along the circum- 
ference of a small circle whose centre is the point of intersection 
in the direction in which is measured. Then it may be shown 
that as we pass from a P curve to a Q curve, P and Q have 



III.] 



OBOMETEICAL ILLUSTRATION. 



89 



opposite signs, and as we pass from a Q curve to a P curve the same 
sign. This is in fact merely Cauchy's rule for the changes of sign 

P 

of -^ when we travel round a radical point. 

23. Let us express the condition that there is no radical point 
on the positive side of the axis of y. This is the geometrical 
illustration of Art. 2. 

Draw a circle of infinite radius, and let it cut the asymptotes 
of the P curve in P^, P^, Pj,...P2^ and the asymptotes of the 
Q curve in Q^, Q^f...Q^» These points alternate with each other. 
Taking only those points which lie^ on the positive side of the 
axis of y, the P and Q curves may be said to begin at these points 
and are to intersect each other only on the negative side of the 
axis of y. The branches of the two curves must therefore remain 
alternate with each other throughout the space on the positive 
side of the axis of y. Their points of intersection with the axis 
of y must be also alternate, and hence if we put x = 0, in the 
equations P=0, Q = 0, and regard them as equations to find y 
the roots of each must separate the roots of the other. 

Conversely, we may show that if the intersections of the two 
branches are alternate on the axis of y, they cannot have intersected 
on that side of the axis of y on which the common intersection of 
all the asymptotes is not. This is the result arrived at in Art. 2. 

24. The following diagrams exhibit the forms of the curve 
p=0, ^ = for a biquadratic. The dotted lines represent the 
asymptotes. 




A biquadratic with four imaginary roots. 
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A biqnadntia with equal imaginai? HMts. 
25. [If a = a: + y V* — 1, we have 

/(2)-P+0/^. 

DifferentiatiDg with respect, firstly to x and secondly to y, we 

It easily foDows that ^ =4^ and ^ = -^ ; so that both 
■^ ace ay - ax ay 

the functions P and § satisfy the equation 



The equation /(a) = gave us t*o curves which we have called 

= Oand Q = 0, In the same way the derived equation /'{«) = 

will give us two other curves, which we may represent by f = 



and Q' = 0. These we may call the derived F and Q curves. 

If as in Art 22 we transform the origin to the point (A, k) we 
have 

P' = CiCosa, + 2(ycoB(5 + oJ+3c/co8(25 + aj) +... 
Q' = Cj sin 0, + 2c,r sin (5 + a^ + 3c,»^ sin (25 + o,) + . . - 
If the origin be at a point of intersection of the curves P = 0, 
P" =0 which is not a double point on the first of these curves, we 
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have Co cos ag = and cos Oi = 0. Hence the tangent to the curve 
P = at the point of intersection is parallel to the axis of x. 
Conversely if we move the origin to any point on the curve P = 
at which the tangent is parallel to the axis of x, we find that the 
curve P = passes through the origin. Hence the derived P 
curve passes through all those points on the curve P = at which 
the tangent is parallel to the axis of x, and all those points on the 
curve § = at which the tangent is parallel to the axis of y. In 
the same way the derived Q curve passes through those points on 
the curve Q = at which the tangent is parallel to the axis of a?, 
and -those points on the curve P=0 at which the tangent is 
parallel to the axis of y. 

If CoCosa^= 0, and Cj =0, the origin is a double point on the 
curve P=0 and the origin also lies on both the curves P' = 0, 
C = 0. So that both the derived curves cut the curves P = and 
Q = in their double pointa In other words, these double points 
are radical points for the derived equation /'(«) =0. In the same 
way, any multiple point on either curve is a multiple point of one 
degree less multiplicity on both the derived curves. 

If a finite straight line AB be drawn parallel to the axis of x 
joining two points A, B on the same or on different branches of 
the curve P = 0, this finite straight line must cut one or more 
branches of the derived P curve. For it is clear that if P vanishes 

at A and P, P which is equal to -j- cannot keep one sign between 

A and P, and must therefore vanish somewhere between A and B. 
If A and B be adjacent points, i.e. if there be no other pomts 
between A and B belonging to the curve P = 0, then the straight 
line AB must cut an odd number of branches of the derived P 
curve. In the same way if a straight line CD be drawn parallel 
to the axis of y joining two adjacent points on a derived Q curve, 
this straight line must also cut an odd number of branches of a 
derived P curve between G and P. 

By considering a tangent as the limit of secant, it again follows 
that if a tangent be drawn to the curve P= parallel to the axis 
of X, the derived P curve must pass through the point of contact. 

Let P be a radical point on the derived curve /'(«) = 0, and let 
it not be a double point on either of the curves P= 0, Q = 0. Let 
a straight line be drawn from R in any direction cutting the 
branches of either of the curves P= 0, ^ = in the points A^ A^ 
&c. Then we may show that 

so that the polar line of R with regard to either of the curves 
P = 0, Q = is at infinity. 
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The positions of the radical points of the derived equation 

f{z)=^0 relatively to any branch or branches of the curves P = 0, 

Q = may be found by the use of Cauchy's theorem. If the 

point spoken of in Art. 1 travel along a branch of the curve P=0, 

it is easy to see that 77 = S^ • If i** travel along a branch of the 

curve ^ = we have ?y = — t" • ^ then any contour be partly 

bounded by branches of these^cnrves, the simplest inspection of 
the points at which the tangents are parallel to the axis will 

determine the changes of sign of y^ as it passes through zero. 

If another part of the contour be an arc of a circle of infinite radius 
whose centre is the origin, the changes of sign through zero will 
be from + to — and their number will be indicated by the number 
of asymptotes of the derived P curve which cut the arc] 

26. The use of Watt's Governor in the steam engine is too 
well known to need description. It has however, as commonly 
used, a great defect. It is sometimes of importance that the 
engine should continue to work at the same rate notwithstanding 
great changes in the resistances. Suppose the load suddenly 
diminished, the engine works quicker, the balls diverging cut off 
the steam, and the engine, after a time, again works uniformly, 
but at a different rate from before. The balls as they open out 
or close in are usually made to describe circles. Let them now be 
coilfetrained to describe some other curve which we may afterwards 
choose so as to correct the above defect. If this curve be a parabola 
and the balls be treated as particles, it is clear from very elemen- 
tary considerations that these will be in relative equilibrium only 
when the engine works at a given rate. This principle is due 
to Huyghens, see Astronomical Notices^ December, 1875. It is 
now proposed to determine the condition of stable oscillation about 
a state of steady motion. 

Two equal rods AB, AB are attached at A by hinges to a 
small ring which can slide smoothly along a vertical axis. The 
ring is attached by a rod to the valve and can thus govern the 
amount of steam admitted. Two equal balls are attached at B 
and By and the centre of gravity O of the rod AB and the ball B 
is constrained to describe some curve. To represent the inertia of 
the engine we shall suppose a horizontal fly-wheel attached to the 
vertical axis whose moment of inertia about the axis is /. Let the 
excess of the action of the steam over the resistance of the load 
be represented by some couple whose moment about the vertical 
axis is f{9), where is the inclination of the rod AO %o the 
vertical, and / is a function which depends on the construction of 



III.] watt's governor. 43 

the engine. Since the steam is cut off when the balls open out, 
it is clear that -^4k^ is negative. 

There may be also some resistances which vary with the 
velocity. Let these be represented by a couple B -^ tending to 
retard the motion round the vertical axis and a couple round A in 
the plane BAS equal to mO -r- , where ^ is the angle the 

vertical plane BAB' makes with a fixed vertical plane. 

Let m be the mass of either sphere and rod; k the radius of 
gyration about an axis through perpendicular to the rod, and k' 
ihs^ about the rod^let l=^AQ» Then thd equation of angular 
momentum gives 

Let the steady motion be given by ^ = a, -^ = w, and let the 
oscillations be represented by ^ = a + a?, ;^ = ** + y« We have then 

/(a) = Bn. 
The equation then reduces to 

{I + 2m (jeTf sin« a + P cos* a)} J + 2m (P + A* - F) sin 2an ~ 

=f(a)x-By. 
This equation may be briefly written 

Let z be the altitude of above some fixed horizontal plane. 
Then if T be the semi vis viva 

2 r= J^'* + 2m {(jfc* + sin»^ + &'* cos*^} ^^ + 2mi^e^ 

If U be the force function, omitting the couple of resistance, 
we have U= — 2mgz, The virtual moment of the couples of 
resistance being 2mCffB0, the Lagrangian equation of motion 

becomes 

d dT dT dU ^ ^^ 

dt dff dd dO 
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Substituting for T and TJ we have when the motion is steady 

!L'(P+jfc«-ifc'«)8in5cos^ = J, 

...|!(j» + ;fc«-ifc'«)sin"^ = «. 

Since Z sin ^ is the distance of Q from the vertical axis, we see 
that the path of G must be a parabola. The semi latus rectum is 

which, we notice, is independent of the radius of the balls. The 
length of this latus rectum must of course be adjusted to suit the 
particular rate at which the engine is intended to work. 

When the system is oscillating about the state of steady motion 
we have, putting 

a« = P + ^(Z" + it»-Jfc'7sin«a, 

and rejecting the squares of x and y, ' 

(ifc» + a»cos"a)^ + (7^-(P + it»-ife'»)sin2any = 0. 

The term a?, it will be noticed, has disappeared from the 
equation. This equation may be briefly written in the form 

Eliminating y from the equations'of motion we have 
^F J + (^(7+ 55^) g + (B(7+ ^i)^ + ZLa? = 0. 

The coefficients are all positive, the necessary and sufficient 
condition of stability is therefore 

{AG^BIC){BG'VEL)>AFEL. 

In some clocks to which Watt's Governor is applied, there is a 
special arrangement which causes G to be much greater than B. 
See the Astronomical Notes, xi., 1851, and the Memoirs of the 
Astronomical Society, Vol. xx. Neglecting therefore B, we have 

GE>FH. 

. :. 2 C^n (P + A;* - Jfc'^ sin 2a >-P(ifc' + a' cos* a). 
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Formation of the equationa of steady motion and of smaU oscillation 
where Lagrang^s method mmf he used. Arts. 1 — 5. 

The equations being all linear the conditions of stability are expressed 
by the character of the roots of a determinantal equation of an even order. 
Art. 6. 

Mode qfeospanding the determinant. Art. 7. 

A msthod of finding the proper co-ordinates to make the coefficients of 
the Lagrcmgianftmction constant. Arts, 8 — 10. 

How the Harmonic oscillations about steady motion differ from those 
about a position of equilibrium. The forces which cause the difference a/re 
of the ' nature of centrifugal forces produced by an imaginary rotation 
about a fixed straight line. Arts, 11 — 19. 

Reduction of the funda/mental determinant to one of fewer rows by the 
elimination of all co-ordinates which do not appear except as differential 
coefficients in the Lagrangian function; unth an example. Arts, 20 — 23. 

Formation of the equations of Motion a/nd of the determinant when 
the geometrical equations contain differential coefficients, so that La/grang^s 
m^ethod cannot be used ; toith an example. Arts, 24 — 27. 

1. Let the system be referred to any co-ordinates f , rj, f, &c. 
The general expression for the kinetic energy is 

where P,. Q, &c. are known functions of f, 17, f, &c. and accents 
have their usual meaning. Let us suppose the system to have 
some motion represented by 

f=/(0. V = F{t), &c. 

and when disturbed, we wish to find the oscillations about this 
motion. To effect this, we put 

?=/(0 + ^> V = F(t)-\-<f>, &c. 
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where 6, <f>, &c. are all small quantities. Substituting and expand- 
ing T in powers of 0, 0, &c., we find 

+ 5/ + 5,f + ... 

+ ^{A^ff' + 2AJ<f>+...) 
+ i^{B^ff' + 2B^ff<f>\+..,) 
+ Cje' + CJ<l>'+CJ'<f> + .,. 

+ &C. 

In the same way we may make an expansion for the Potential 
Energy of the forces, viz, 

+ ^{E,,e' + 2EJ<f> + ...), 

when these two functions are given the whole dynamical system 
and the forces are known; and we may form the equations of 
motion by Lagrange's method. 

2. We shall here however limit the question by supposing 
that the motion about which the system is oscillating is what has 
been called in Chap. i. steady. The analytical peculiarity of such 
a motion is that when referred to proper co-ordinates, every coefiB- 
cient in each of these two series is constant, i.e. independent of t 
As already explained the physical peculiarities are that the vis 
viva is constant throughout the steady motion and the same oscil* 
lations follow from the same disturbance at whatever instant it 
may be applied to the motion. A method of discovering the 
proper co-ordinates, if unknown, wiU be given a little further on. 

3. In order to form the equations of motion we must now 
substitute in Lagrange's equations 

. d_dT_dT dV_ 
dtdff d9'^~d0'~ ' 

&c. = 0, 

rejecting all the squares of small quantities. The steady motion 
being given by 0, 0, &c. all zero, each of these must be satisfied 
when we omit the terms containing 0, ^, &c. We thus obtain the 
equations of steady motion, viz. 

A = K 

&C. = &C. 
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These equations may be simply formed in any case by the 
following rule. Putting L=T^V, so that L is the diflference 
between the kinetic and potential energies, expand L in powers of 
the co-ordinates 6, <}>, &C; regarding ff, ^', &c. as zero. The required 
relations are obtained by equating the coefficients of the first 
powers to zero. This rule may be also expressed thus. Let L be 
the general expression for the excess of the kinetic energy over the 
potential energy of a dynamical system in terms of its n co-ordinates 
f, rj, &c. Let this system be moving in steady motion with constant 

values of -^ y-TL, &c.* Then, svbstituting these constant values in 

the general eocpression for L, the relations between the constants of 
steady motion are given by 

dL dL ^ « 

In this way we obtain in general as many equations as there 
are co-ordinates. Usually the coefficients in the expression for T 
are constant because some of the co-ordinates are constant in the 
state of steady motion, and the other co-ordinates appear in the 
expressions for T and V only as differential coefficients. In such 
cases we have clearly fewer equations than co-ordinates to con- 
nect the constants of steady motion. We have then a system of 
possible steady motions which we may conveniently term parallel 
steady motions. 

4. To obtain the' equations to the oscillatory motion, we retain 
the first powers of 0, </>, &c. We thus obtain a series of equations 
of which the following is a specimen : 

i + =0. 

To solve these we write 

0=-M/^, (l> = M,e'^, ylr = M/^, &c. 

Substituting we obtain on eliminating the ratios M^ : M^ : M^ &c. 
a determinantal equation, viz. 

* Since we may in Art. 1 change the co-ordinates from ^, 17, f, &c. to ^j, tf^, f j, 
Ac, where ^=/(^j) 7i = F(rfj) Ac, it is clear that the steady motion can be always 
expressed by constant values of the differential coefficients of the co-ordinates. 
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B^m' — Aji + Ejj^ 


B„m*-A^ + E^ 


Bam'-A„ + E„ 

+ (Cn - Cu) TO 


&c. 




B^ -Aa + E^ 


B„m*-A„ + E„ 


&c. 


B„m'-A^+E,, 


B^m'-A^ + E^ 


B„m'-A^ + E^ 


&c. 


&c. 


&c 


&c. 


&c. 



= 0. 



This equation will be referred to as the Determinantal equa- 
tion, 

5, I£ we refer to the equation formed by this determinant 
and read it in horizontal lines, we have of course the several 
equations of motion, each term being the coefficient of 0, <f>, yfr, &c. 
in order. In this form the equations may be reproduced by the 
following easy rule. 

Taking the expression for T— V as given in Art. 1, let us 
consider only the terms of the second order, those of the first 
having been already used to determine the steady motion as ex- 
plained in Art. 3. Separate from the rest, the even powers of 

0> &> <l>, <f>\ &c. and write for ff\ 0'<f>', &c. - L^O", - jD'^^, &c., so 

ft 
that D will stand either for -r- or for the m in the determinant 

at 

of Art. 4, when we write = M^, <f> = M^eT*, &c. Let the sum 
of these terms be called P, so that 

P = i {A,, - E,, - B,,L^ ^ + (J., - E,, - B,,I>') e<l> + &c. 

Let the remaining portion of the terms of T— F, viz. those con- 
taining both 0, (f>, &c. and ff, <f>, &c., be called Q, so that 

Then the several equations may be formed from the rule 

dff^ d0 ^ dff~^ 

d(f) d<j> d<fl 

&c. = J 

In applying this rule no accented letters will occur except in 
the second term of each equation. If we wish D to stand for m 
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in the determinant, we must regard 0^ ff>', &c. as abbreviations for 
1)0, D(f>, &c. If we wish to use the equations themselves, we 
replace 1)0, D^, &c. by 0\ <^', &c, 

fThe determinant may also be found by another rule. Taking 
aslLfore only the termsV th. second oLr in the Lagrangiaf 
function i/ = 2^— F, let us separate the terms of- the form 

Q = CJff + CJfi; + GJ^4> + &c. 

In the remaining terms put ^ = ^m V— 1, 0' = ^m V— 1, and 
so on, and write down the discriminant If the system oscillates 
about a position of equilibrium the terms represented by Q are 
absent and the discriminant thus formed will be the determinantal 
equation giving the required values of m. But if the system 
oscillate about a state of steady motion we must modify the dis- 
criminant by adding some quantity derived from Q to each term. 
To find this, write above the columns ff, (f>, &c. and before the 
rows 0, <^, &c. Consider any term, say the term in the column 
<f)' and row 0. We must add to that term {G^^— C^^ m, where 
(7,a — C^j is the excess of the coefficient of 4/0 above the coefficient 
of <f>ff in the expression for Q, Since the determinant is un- 
changed by writing — m for m, we may, if preferred, add the 
excess of the coefficient of ff<f> above the coefficient of 0<f} provided 
we adhere to one order throughout.] 

6. K in the determinantal equation we write — m for m, the 
rows of the new detenninant are the same as the columns of the 
old, so that the determinant is unaltered. When expanded, we 
shall have an equation which contains only even powers of m. 

The condition of dynamical stability is that the roots of this 
equation should be all of the form 

If we write m^ = —p\ the roots of the transformed equation 
must be aH real. 

In the case of equal roots of the form m= + 0V— 1, it has 
been shown in Art. 5 of Chap. I. that it is necessary for stability 
that the proper number of minors in this determinant should 
vanish. If there be two equal roots, these roots must make all 
the first minors zero ; if three equal roots, all the first and second 
minors must vanish, and so on. 

7. When the system depends on many co-ordinates the labour 
of expanding this determinant is often considerable. Methods of 
evading this in certain cases will be given in the following chap- 

K. A. 4 
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[CHAP. 



ters. But when the development is necessary we may proceed in 
the following manner. Let the determinant be written 



-B,,m*+an 



&c. 



&c. 



+ F^m 

&c. 



and let 



0" .^ 

Pi flP* 

We know that the determinant when expanded is of an even 
order, hence all odd powers of m must finally vanish. Let us 
expand the determinant in powers of the J"s. The first term is 
the discriminant of Bm^ + a, this term is independent of the -Ps. 
The terms which contain the first powers of the F's are obtained 
by erasing any one line of this discriminant and replacing it by 
the corresponaing F terms. But these terms all vanish and we 
need not describe them minutely. The terms containing the 
products and squares of the F's may be obtained by erasing every 
two lines of the discriminant and replacing them by the corre- 
sponding J* terms. Thus if we erase the two first lines we have 
the determinant 




-Fjn 



IS' 



&c. 



F^^m 


&c. 



F^m &c. 
Fjn &c. 
Bjjn^ + a^ &c. 
&c. 



and so on for all the other rows taken two and two. The terms 
which contain the cubes and all odd powers of the ^s vanish, 
while the terms which contain the fourth powers may be obtained 
by erasing four lines of the discriminant and replacing them by 
the corresponding Fs. 

When the determinant has been expanded, we have an equa- 
tion of an even order to find the values of m. We may therefore 
employ the short method of Art. 5, Chap, iii., to obtain the 
Sturmian functions. 

8. Necessary and sufficient tests of the stability of the motion 
of a system of bodies are given in the preceding pages. But it 
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• 

is assumed, as explained in Art. 2, that the co-ordinates have been 
properly chosen. They are supposed to have been so chosen that 
the coefficients in the expanded Lagrangian function are all 
constants. When this is not the case we must discover the 
proper co-ordinates to which the system must be referred before 
we can apply the test of stability. But when the motion is steady 
this is not difficult. 

There are obviously many such systems of co-ordinates, and 
one set may generally be found by a simple examination of the 
steady motion. If there are any quantities which are constant 
during the steady motion, they will often serve for some of the 
co-ordinates. . Others may be found by considering what quan- 
tities appear only as diflferential coefficients or velocities. Practi- 
cally these will be the most convenient methods of discovering 
proper co-ordinates, since no further change will then be necessary 
and we may at once form the determinant of stability. But if 
these methods fail we may adopt the following analytical method 
of transforming (where possible) the general Lagrangian function 
with variable coefficients into one with constant coefficients. 

• 

9. Let the Lagrangian function be 

'i-^A^^0' + A^j9(f>+... 

+ ^B,,ff' + B,,ff<l>' + ... - 

where the coefficients are all functions of t and the co-ordinates 
0, 4>y &c« have been so chosen as to vanish along the steady motion. 
We have therefore for the steady motion 

The oscillations about the steady motion are given by the 
terms of the second order. Our present object is to transform 
these to others with constant coefficients by the following substi- 
tutions : 

6 ^p^x +PJJ -{-p^z + &c., 

&c. = &c., 
where the p's, j's, &c. are functions of the time at our disposal 

Substituting and equating the coefficients of x^, y'*, &c, to 
unity, we have as many equations of the form 

i5,^« + i5„2' + 5„pj+... = l (1) 

4—2 
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as there are co-ordinates. Equating the coefficients of the pro- 
ducts xy ^ x*z\ &c. to zero we get » — 5— equations of the form 

Ail>ii^. + ^MM8 +-S,, (m«+JP2?i) + -•• = (2)> 

supposing that there are n co-ordinates. 

Equating to the constants a^, a,, ... the coefficients oixx\ yy', 
&c. having subtracted the diflferential coefficients of (1) we have 
n equations of the form 

(C^«-iA.')i''+(Ci.-i^a')2'+(^«+C'.i-A,')2'2+-=«-(3)- 
Adding the coefficients of xy' and xy and subtracting the 

differential coefficients of (2) we have n — ^ — equations of the 

form 

Equations (1), (2) and (4) give n? equations to find the t? 
quantities p^p^^ &c., q^q^, &c. The solution of these equations is a 
purely geometrical problem. If we construct the two quadrics 

\B,,0'-\-\BJ?-]rBJ<t> + ... = 1, 

(On-i5'n)^ + (a,-iB'J^«-f((7«+(7.,--B'J^</>-l-...=l, 
and refer them to their common conjugate diameters, by writing 

e=p^x+pj/+... 
<t>-qiX + qj)-\- ... 
&c. = &c., 

making the first quadric to become what we may call a sphere by 
projection; the values oip^p^ &c., q^q^ &c. thus found are the values 
required to make some of the coefficients in the Lagrangian 
function become constant. These values must of course make all 
the other coefficients of the second order m the Lagrangian function 
constants also, and thus we haver^(n + ]) analytical conditions 
that the motion should be steady. 

It might be supposed that greater generality would be obtained 
lay replacing the zero's of equations (2) and (4) or the unities of 
(1) by arbitrary constants. This may be convenient in practice, 
but as we know that by a subsequent real change with constant 
values of p^p^y &c., q^q^, &c., we can render them zero or unity, it 
simplifies the argument to perform the two transformations at 
once. 
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10, The geometrical problem just alluded to admits of a real 
solution whenever one quadric can be projected by a real projection 
into a sphere. The problem then becomes that of fiDding the 
principal axes of the other. This is just our case, since the ex* 
pression 

is necessarily positive for all values of ff^ . 

It is unnecessary to describe here the mode of solving this 
problem. It is sufficient to say that it may be reduced to the 
solution of the symmetrical determinantal equation 

= 0, 



5ii — \2?ji. 


Jt/]g — A,X/22 • • • 


B^ — ^-Dj2» 


X/22 ~" '''•^28 • • • 







where Du, 2)^, &c. are the coefficients of ^, ^^, &c. in the second 
of the quadrics. The roots of this equation are known to be real 
when the suppositions just mentioned are satisfied. 

11. In order to examine the fundamental determinant in 
Art 4 a little more closely, let us suppose it reduced to depend on 
three co-ordinates. We may then have the advantage of a geo- 
metrical analogy. Let the co-ordinates be f ^, 5> ai3<J kt the 
equations of motion be written 



(i.4-»)f-r"5,"'Vf": 



de ^" 



. F— 
dt ' ^ dt 





d* . . „ <P 



( _4 y\ ,^| j'^(*-^-^")'=»- 

Let a geometrical point P move in space so that its co-or- 
dinates referred to any axes are ^, rj, ^. Then the position and 
motion of the point exactly give us the position and motion of 
the system. 

12. Looking at the equations of motion just written down we 
see that they are similar to those which give the oscillations of a 
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system about a position of equilibrium, but that there are in 

addition terms E -^, JF-f , &c. The general effect of these terms, 

at at 

as will appear from what follows in the subsequent chapter, is to 

increase the stability. If we transpose these terms to the other 

sides of the equations, we may regard them as impressed forces 

acting on the system, whose resolved parts in the directions of the 

axes f , rj, f, are 

at at 



^=4^-4'. 



We see at once that 



EX+FY+OZ^O 





dt dt dt 



so that these forces are at once orthogonal to the path of the 
representative point P and also orthogonal to the straight line 
whose direction cosines are proportional to E^ F, 0, These forces 
are therefore of the nature of centrifugal forces^ as if they were 
produced by the rotation of the system about this straight line. 

13. We may show that the straight line {E, jP, Q) is fixed in 
space. To prove this, let us transform our co-ordinates from f , 17, f 
to a?, y, 2r, where a?, y, z are connected with f , 1;, f by any linear 
relations, such as 

f = a^a? + 1 J/ + c^z\ 
rj^a^-^rhjf+c^ 

Let the portion of the Lagrangian function under consideration 
(Art. 1) be 

then G = C,, ~ (7,„ E^ C^ ^ C^, F^G,,^ €,,. 

Substituting for f , ^ ? their values in terms of x, y, z, we find that 
the diflference between the coefficients of xy' and xy is 

G'=0 



«.«« 


+ E 


«.a. 


+ F 


»a<*i 


hh 




Kh 




b,h\ 



♦: 



I 
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with similar equations for E and F. But if /i be the determi- 
nant of transformation 



«»«« 


+ ? 


«.«. 


+ '7 


»•«! 


\K 




KK 




J. Jl 



with similar equations for x and y. The ra^io5 of ^, J?) (7 are 
therefore transformed as if they were co-ordinates. If the trans- 
formation be a real transformation of Cartesian co-ordinates, let 
lengths each equal to unity be measured from the origin along the 
axes Of, Oi;, 0£ thus forming a tetrahedron whose volume is F. 
Let a similar construction be made for the new axes, forming a 

tetrahedron of volume Y. Then* fi — 'Tr* Hence the quantities 

E F G 

pr, Tf> 'n- ^^7 be transformed as if they were lengths measured 

E F Q-' 
along the axes and become y/, y?, y? • If both systems of co- 
ordinates are rectangular we have F= x , T^ = q ♦ 

jr w Q- 

Let ft) he the resultant of y^^ ^ , y , then ft) may he regarded 

as a fixed length measured from the origin along a straight line 
fixed in space. 

Let V be the velocity of the representative particle, the angle 
between the direction of this velocity and the axis whose direction 
cosines are proportional to B, F, O. Then the resultant of the 
forces Xy F, Z is easily seen to be 2t; Fo) sin acting perpendicular 
to the axis and to the direction of the motion. We might call the 
straight line (JEFG) the axis of the centrifugal forces. 

[* Let (^1 i7i ^i), (f J i7j fj), (fg i/g fj) be the oo-ordinates of three points A, -B, C re- 
ferred to any oblique co-ordinates. Let ns find the yolume V of the tetrahedron 
of which these and the origin are the angular points. Since the yolume vanishes 
when any angular point as C Ues in the plane containing the origin and the other 
two A, Bf the expression for the volume must contain the factor 



fi= 



h ^t U 

Vi Vt Vz 



The volume is evidently an integral rational function of the co-ordinates when 
the axes are rectangular and the plane ^0^ is taken as the plane of xy, it 
easily follows that this is true for all axes. Since this function cannot be of more 
than the third order, we have V'=Mfjk, where M is independent of the co-ordinates 
ot A, B, C When the points A, B, C are on the axes at unit distances from the 
origin, let V be the volume of the tetrahedron. In this case /t=l, and .*. M=V, 
We have therefore in all cases V = 7/*. 

In the text, let the extremities of the unit lengths measured along the axes of 
«, y, z be called A, B, C. Then tiie (^17^) co-ordinates of ^, B, C7 are (a^, a^ a,), 
(6 J, &,, &,), (Cp c,, c,), respectively. Hence by what has just been said K'=K/t.] 
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14. The expressions for the co-ordinates in terms of the time 
wiU in general contain as many periodic fonctions as there are 
co-ordinates. If the initial conditions are such that each contains 
one and the same periodic function, the motion recurs after a con- 
stant interval and the system is said to be performing a dmple or 
harmonic oscillation. 

If the system be oscillating about a position of equilibrium, 
with a Lagrangian function 

we know* that the hartnonic oscillations are represented by rectili- 
near motions of the representative particle, and that these are 
along the common conjugate diameters of the two quadrics ^ ' 
equations are 

where a and B are two constants chosen to make the quadrics real. 
Let us consider what are the harmonic paths of the representative 
point when the system is oscillating about a state of steady motion. 

In any harmonic vibration we have f = i.co8 (Xi + a) with 
similar equations for rj and (f. Hence 

df--^ir ^ A'^, -^ A,f. 

Substitute these in the equations of Art. 11. Differentiate and 
substitute again. Multiply by f, 17, 5", add and integrate, we obtain 

where c is some constant. The harmonic path lies on this quadric, 
which has a common set of conjugate diameters with the two 
quadrics a and 6. 

If we resume the result of the substitution of j5 > &c. in the 

dv 

equations of Art. 11, and multiply by -F, Fy respectively and add, 

we obtain 

1{B^E+ Bi^+B^G) i + &c.] \» + [{A^E+ ...) f + &c.] = 0, 

* A short paragraph in Thomson and Tait's Natural Philosophy, page 273, is 
the only notice of this which the author has discoyered. 
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which is a plane, and is diametral to the straight line [EFG) with 
regard to the quadric c. 

The harmonic paths are therefore ellipses. The three harmonic 
planes are diametral to the same straight line and this straight 
line is fixed in space, being the axis of the centrifugal forces. 

If we eliminate \ between the equations to the plane and the 
quadric c, we get a cubic surface on which the three harmonic 
conies lie. 

If Ey F, are zero^ which is the case when the system oscillates 
about a position of equilibrium, the quadric a becomes a cylinder. 
This may be conveniently shown by referring the system to such 
co-ordinates that the coefficients B^^ B^, B^, A-^, A^^ A„ are all 
zero. In this case the diametral plane of every straight line passes 
through the axis of the cylinder. The harmonic oscillations are 
^ therefore rectilinear. 

If B be the length of that semidiameter of the quadric (c) 
which is parallel to the fixed straight line (E, F, G), it may be 
shown that the 

Product of the axes ] R 2c 

of the quadric c ) JW^TF^^^^ * ^ * 

If ^, F, O are all zero, and their ratio is indeterminate, B is 
any diameter. Hence one of the axes of the quadric (c) must be 
infinite and the quadric will be a cylinder. 

[If the quadric (c) be a cylinder and E^F^ O are not all zero, 
we must have either \ zero or B infinite. In the latter case the 
axis of the cylinder will coincide with the axis of the centrifugal 
forces.] 

The quadric (c) has also the following geometrical property, 
Xet the lengths of semidiameters of the quadrics (a) and (6) drawn 
parallel to the axis of the centrifugal forces be p and p . Through 
the intersection of'these quadrics describe a quadric so that the 



Product of) 2VaJ 

its axes 

'"" ' ■ ■ '^ A/ I 

This quadric is similar to the quadric (c). 



(7*-?)^ 



V V iPAi? pV 



15. The introduction of the representative point to exhi1;>it 
the motion of a system may appear somewhat artificial. If how- 
ever we properly choose the co-ordinates the particle moves exactly 
as a free particle, and we might reduce the problem of finding the 
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oscillations of a system to a problem in Dynamics of a particle. 
Refer the quadric 

to its principal axes and let the equation thus changed be 

B^ ^+ J5j,' — + ... = 6. 

Since jBu , B^, &c. are positive quantities, we may put 

^B'^ix-^x, V:^'^, = y, &C. 

The quadric has thus been "projected" into a sphere. Let a;, y, z 
be now chosen as the co-ordinates of the system and let the 
Lagrangian function be expressed in the form 

"t" V/jjfl/iC "J" • • • 

the terms of the first degree being omitted as not necessary to our 
present purpose. The three equations of motion at the beginning 
of Art. 11 take the form 

&c. = 0, 

which are the three equations of motion of a free single particle 
whose co-ordinates are a?, y, z under the action of forces whose force 
function is 

and a force acting perpendicular to the path and also perpendicu- 
lar to a fixed straight line, the force being proportional to the 
velocity. 

16. As an illustration of this theory, let us here make a short 
digression. However the particles of light may oscillate, whether 
in a rotatory or linear manner, we know the motion is related to a 
certain plane called the plane of polarization. It may be shown 
that any harmonic oscillation about a position of equilibrium may 
be represented by a rectilinear oscillation of the representative 
particle. Let us represent the motion at any point of the ether 
by a rectilinear oscillation in a direction perpendicular to the 
plane of polarization. This would be FresneFs Vibration. The 
representative particle, as just shown, would not necessarily move 
as if it were a free single particle. But let us assume (and a proof 
is not necessary to our present purpose) that when the oscillation 
is drawn as above described the motion in the plane of the front 
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is the same as that of a free particle, while that perpendicular is 
not free. On this assumption we see that Fresnel in his theory 
of double refraction is justified in taking actual instead of relative 
displacements, for it is the representative particle he is considering. 
He also neglects the force normal to the front, for the particle 
moves as a free particle only in the plane of the front. These 
general remarks are not meant to explain Fresnel's theory, but 
merely to show how the representative particle may be used to 
replace a complicated motion. 

17. [When a system is performing a harmonic oscillation 
about a state of steady motion or about a position of equilibrium, 
the motion repeats itself continually at a constant period, that is 
to say, the values of the co-ordinates recur at this interval. This 
is the chief peculiarity of a harmonic oscillation. When the 
oscillation is about a position of equilibrium, the representative 
particle oscillates in a straight line whose middle point repre- 
sents the position of equilibrium. Thus the system passes through 
the position of equilibrium twice in each complete oscillation. 
When the oscillation is about a state of steady motion the path 
of the representative particle is an ellipse whose centre is at the 
point occupied by the system in steady motion at the same in- 
stant. Thus the system does not in general ever coincide with 
the simultaneous position of the system in the undisturbed or 
steady motion. When a system is disturbed by a small impulse 
from a state of steady motion, it will in general describe a. com- 
pound oscillation made up of at least two harmonic oscillations, 
at the instant of disturbance these two neutralize each other so 
that in the disturbed and undisturbed motions two simultaneous 
positions are coincident. But it is clear that this cannot occur 
again unless either the periods of the two harmonics are com- 
mensurable or the period of one of them is infinite.] 

18. [In some cases the ellipse degenerates into a straight line. 
Thus if the quadric (c) be a cylinder the diametral plane of the 
axis of the centrifugal forces will pass through the axis of the 
cylinder, and thus the harmonic oscillation corresponding to this 
particular value of \ will be rectilinear. In this case the system 
twice in each oscillation passes through the position it would have 
occupied at the same instant in the undisturbed motion. 

The quadric (c) has a common set of conjugate diameters with 
the quadrics (a) and (J). Hence if (c) be a cylinder, its axis must 
be parallel to one of the three common conjugate diameters of (a) 
and (b). If we refer the quadrics (a) and (6) to their common 
conjugate diameters, they take the form 



7 ^ 
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The cylinder which passes through their intersection and has 
its axis parallel to the diameter f is found by eliminating f* between 
these equations. We see therefore that -Bgg'X* + A J = 0. If then 
the axis of the cylinder cut the quadrics (a) and (6) in D and iX 
respectively, we find that for this oscillation 

, on^ -a 

~ OD^ ' b ' 

It has already been shown that when this value is finite, the 
direction of ODU is along the axis of the centrifugal forces.] 

19. [In some cases two or more of the values of X are zero. 
In these cases the co-ordinates will have terms of the form nt + e, 
where n and e are two small constants. When, as explained in 
Art. 3 of this Chapter, there are several parallsl states of steady 
motion, these terms imply that the motion is stable about a state 
of steady motion very nearly the same as the undisturbed motion 
but not coincident with it. The actual undisturbed motion, unless 
n is zero, is unstable in the sense that if a proper disturbance 
be given to the system, the system will depart widely from the 
positions it would have simultaneously occupied in the undisturbed 
motion.] 

20. In many cases of small oscillations it will be found that 
the Lagrangian function T— F is not a function of some of the 
co-ordinates as 6, <f>, &c. though it is a function of their diflferential 
coefficients ffy <f>', &c. In such cases the steady motion will be 
usually given by constant values of these differential coefficients, 
while the other co-ordinates as f, rj, &c. are also constant. It is 
evident that the determinantal equation of Art. 4 is needlessly 
complicated. It is clear that there will be as many pairs of roots 
equal to zero as there are co-ordinates 0, 0, &c. It will be an 
advantage to eliminate d', 4>, &c. altogether from the Lagrangian 
function, and to find the remaining roots by operating only with 
the co-ordinates ^, rj, &c. We shall thus obtain a determinant with 
just as many rows as there are co-ordinates of the kind f , rj, &c. 

Let ij be the Lagrangian function expressed as a function of 
ff<j>\ &c. ^, ^V> &c. Let L^ be its value when 0'<t/ are elimi- 
nated, so that i/j is a function of ^, ^i]', &c. only. To effect this 
elimination we have the integrals 

dT_ dT _ . 

where Cj, c,, &c. are constants. Then 

dL^^dL, dL, dff dL, d£ . 
If' dJ^^ dff 'd^ '^d^' rfr 
__ dij dff d<f> 
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dL,_dL dff_ d^ 

d^~ d^ ^^d^'^'^^d^'^'" 

•'• dt df d^~dt rff ' d? ■•" *'» \dt df d^J "*" **^- 
But ^^^^-0 

Hence if we take 

jrn jm 

and eliminate 0, <f> hy help of the integrals j^ = ^xy ^'~^8» 

we may treat 1/ just as we do the Lagrangian function L. The 
equations giving the small oscillations about the steady motion 
wUl he 

d dL' dL'_ 

dtdr""d|""""' '^^•■"^• 

The function L' may be called the modified Lagrangian fwnction, 

[It should be noticed that this is equivalent to a partial use of 
Hamilton's transformation of Lagrange's equations. Sir W. R. 
Hamilton eliminates all the differential coefficients ffy 0', &c. by 

the help of equations of the form jff'='^> ;t . , = v, &c. where w, v, 
&c. are made to be new variables*. In our transformation only 

* The Hamiltonian transformation of Lagrange's equations bears a remarkable 
analogy to the transformation of Beciprocation in Geometry. This may be shown 
in the following manner. 

When the system has three co-ordinates 6^ 0, ^, We may regard ^', 4/, ^ as the 
Cartesian co-ordinates of a representative point P. The position and path of P wiU 
exhibit to the eye and will determine the motion of the system. Let u, v, \d be the 
Hamiltonian yariables, so that 

dr, dT. AT. 

«=dtf" "=5^/' ""d^' 

where Tj is the semi vis viva expressed as a function of d, 0, ^, ^, ^^ yj/. Then 
Uy Vf w may be regarded as^the co-ordinates of another point Q whose position and 
path will also determine the motion of the system. 

If the semi vis viva be given by the general expression 

it is dear that the point P always lies on the quadric T^^U where U is the force 
lonction and the co-ordinates 0y 0, \ff have their instantaneous values. The point 
Q must i^erefore lie on another quadric which is the polar reciprocal of the first 
with regard to a sphere whose centre is at the origin and whose radius is equal to 
^2 U, The equation to the reciprocal quadric is therefore 






u V w 

** -^11 -^IS -^18 

w A12 A^ A^ 



u, 
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those new Tariables are introduced which would be constants in 
Sir W. R. Hamilton's transformation. 

This remark suggests an extension of the process. If Z be a 
function of 0y <f), &c. as well as of ff, <l>\ &c. the quantities c^, c^, &c. 
will not be constants. We express this by writing u, v, &c. instead 
of Cj, c,, &c. Suppose we wish to eliminate some of the difierential 
coefficients, viz. a, <l>\ &c. and to retain the remaining ones, viz. 
f , r/, &c. If we put 

we may easily show as in the preceding page that 

d dL dL /\ o /\ 

where A is the determinant, caUed the discriminant, which may be formed from 
the determinant jnst written down by omitting the first row and the first column. 

This is a general expression for the Hamiltonian function and agrees with that 
which may be dednced from the result in Art. 21, when aU the variables are trans- 
formed by the Hamiltonian process. 

Since the polar reciprocal of the polar reciprocal is the original quadric, it 
follows that 

0^-^^^ ,6'-£I» O/-"^^' 
^-dH' "^-1^' ^-^' 

which are three of the six Hamiltonian equations. 

We may also show geometricaUy that if the coefficients of J\ be functions of 

any quantity 6, then -^= - -3q • To prove this we notice that if ie, y, « be the 

co-ordinates of a point P, situated on a radius vector OP* of a quadrio 

(0P\* 
jYp, \ . The quadrics 

Ti = l and 7,=1 may be regarded as polar reciprocals of each other with regard 

to a sphere whose radius is ^ and whose centre is the common centre of the 
two quadrics. Let P be any point on the quadrio 2\=1, and let the radius vector 
be produced to Z so that OP . 0Z=2, then the quadric T,=l touches a plane 
drawn through Z perpendicular to OP and Q is the point of contact. Let these 
quadrics be slightly altered in consequence of a variation of 9, so that their 
equations are now Ti-\'dT^=l and T^+dT^=l, Let OP and OQ produced cut 
these quadrics respectively in P^ and Q\ Then 



T,+dT,=(^^y. r.+dr.=(^y. 



Now if Z' be a point on OP produced so that OP . 0Z'= OP . OZ, the quadric 
T^+dT^'wHl touch the plane drawn through Z' in some point q near Q^. The 
point Q' will therefore lie very nearly in the tangent phme, so that by similar 
triangles 

OQOZ _ OP 

0Q:~0Z'~ OP' 

Since each of these ratios is indefinitely nearly equal to unity, it follows that 
dTi=-dT^. 

If we put L=Ti + U and H=T2- U, Lagrange's equations may be written in 

the forms «'= -j;r , »'=:?-, tt/=-77 . Hence we have 
d$' d4> d^p 

, dH , dH , dH 

-"=5^' -^=d0' -«^=5^' 

which are the remaining three of the Hamiltonian equations. 



M^ 



IV.] • 



LAGRANGIAN FUNCTION. 



C3 



We have thus as many equations of the Lagrangian form as 
there are variables f , rf, &c. Also since w = jg? , &c. we have by 
differentiation 

with similar equations for <l>\ &c. By Lagrange's equations we 



obtain 



de 



= u\ &c. 



Thus we have as many sets of equations of the Hamiltonian 
form as there are variables 0^ <f>, &c.] 

21. We may effect this elimination once for all and find a 
definite expression for L\ 

Let the kinetic energy be 

T = T^ -^ + T$^ ff<f> + . . • 

Then the integrals used will be- 

2^w^ + 2V^^' + ... =Ci— Tijff' — T^rf — ... 

&c. = &c. 

For the sake of brevity let us call the right-hand members of 
these equations Cj — X, c^^-Yj &c. Since T is a homogeneous 
function, we have 

.'. z'=rtt^+!r^r^'+&c.-F 
-i^(ci-z)-if(c,-r)-&c.j 

If we substitute in the second line the values of ff, (f/, &c, found 
by solving the integrals just written down, we have 

L'=TtX+Tt,^r,' + &c.-r 



+ 



2A 



<j,-Z c-F 

Ci — X Tm 2V^ 
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where A is the discritnmant of the terms in T, which contain only 
e, f , &c., and may be derived from the determinant just written 
down by omitting the first row and the first column. 

We may expand this determinant and write it in the form 



1 

+ 2A 





• • • • • ■ 


• • • 

• • • 

• • • 


1 

+ 2A 



X 

y 

• • • 




y 

• ■ • 




1 





X 


y ... 








A 


* • 1 


* • • • 


Tf^ ... 

A^ • • • 
••• ••• 


) 




1 - - 

where X, Y, &c. stand for 














I 
I 

&c 


r 


&c. 


1- 





The first of these three determinants will contain only the 
constants c^, c^ &c., and the co-ordinates ^, 1;, &c. The second will 
not contain Cj, Cj, &c. but will be a quadratic function of ^, 17', &c. 
The last determinant will contain terms of the form |', 17' with 
variable coefficients which may also be functions of c^, Cj... 

Since ff, ryq\ &c. are all small quantities, it is clear that this 
expression for L' when expanded will take a form precisely similar 
to that given in Art. 2, only that we have fewer variables to deal 
with. 

22. As an example, let us consider the following problem. 

A body has a point O which is in one of the principal a>xes at 
the centre of gravity G fixed in space. The body is in steady 
motion rotating with angular velocity n about OG which is vertical. 
Find the conditions that the motion may be stable. 

Let OA, OB, 0(7 be the principal axes at and let OG co- 
incide with the vertical OZ in steady motion. Let f , r} be the 
direction cosines of the vertical OZ referred to OA, OB. Let 
6)1, Wj, 6)3 be the angular velocities about the principal axes at (?. 
Then to the first order 
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Let be the angle ZOC, ^jr the angle the plaae ZO (7 makes 
with a plane ZOX fixed in space and ^ the angle it makes with 
the plane AOG fixed in the body. Then 

f = — sin^cos^' 
17 = sin ^ sin <j>) ' 

Q>3 = 0' + '^' cos 

=f+(x'-^')(i-f). 

putting X = ^ + V^. We easily find 

«.=x-x'^'-i(^'-l''»)- 

K then J, £, (7 be the principal moments of inertia at 0, the 
Lagrangian fimction is 

where M is the mass of the body, and h = 00, 

Since x ^ absent from the equation we have the integral 

dL_ 

which gives 

X =n + terms of second order. 

Hence 

L' ssL^ Cnx 

= const + ^,'« + |r 

+ {{A-C)»* + Mgh}^+[iB-C)n*+Mgh}'^ 

Using this as the Lagrangian function we easily find 

(A - C)n* + Mgh - Bm\ (A + B-C)nm =0. 

-(A+B-C)nm, {B - C)n* + Mgh - Am 

B.A. 
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The roots of this equation to find m must for stability be of 

the form +)9\/ — 1. Putting m" = — V we have a quadratic to 
find \\ The roots of this quadratic must be real and positive. 

li A=B, as in the case of a top spinning with its axis vertical, 
we have 

^ (2A-0 . / G'n'' - 4iAMgh) 

^=±i-2Z-^±\/ 2A-\' 

The motion is stable or unstable according as Cn^ is greater or 
less than i^AMgh, If (7W = 4iAMgh, the equation has equal roots 
and as the first minors are not zero the motion is unstable. 

23. [As another es:ample of the use of the modified Lagrangian 
function, let us consider a case discussed by Prof. Ball in the 
Notices of the Royal Astronomical Society for Marchy 1877. In a. 
problem in Physical Astronomy, we want the relative co-ordinates 
of the system, while its absolute motion in space does not concern 
us. Lagrange's equations involve both the relative and absolute 
co-ordinates, and are therefore not particularly well adapted for 
such problems. By using the modified Lagrangian function, we 
may eliminate the absolute co-ordinates. 

Let the system have n co-ordinates, let us choose as three of 
them the co-ordinates of the centre of gravity of the whole system, 
viz. Of (f>, y^. There will then remain w — 3 co-ordinates which are 
independent of these. Let T' be the kinetic energy of the system 
relative to its centre of gravity, V the potential energy, M the 
whole mass. Then the Lagrangian function is 

In problems in Physical Astronomy the potential energy is a 
function only of the relative positions of the bodies, and is there- 
fore independent of 0, ^, -^ and their diflferential coefficients. We 
have therefore 

dL dL ^ dL __ 

dff'^^^' df"^«' rff"^»' 

Hence the modified Lagrangian function is 

i' = y — r— a constant. 

It is independently clear that we might take this as the 
Lagrangian function, for the first three terms of L do not enter into 
any one of the Lagrangian equations, except the three formed by 
differentiating with regard to ff, <f}\ yft'. 

The function T" is made up of two parts, (1) the kinetic energies 
of the rotations of the bodies about their centres of gravity, ^hich 
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we may call 7/, and (2) the relative kinetic energies of the several 
bodies, each collected at its centre of gravity, -which we may call 
Tg'. Let m^, m^, &c. be these masses; x^, x^, &c. the abscissdB of 
their centres of gravity referred to the centre of gravity of the 
whole as origin. Then^ accents denoting differential coefficients 
with regard to the time, we have 

m^a^ + m^x^ + &c. = 0. 

Let us square this and write 

2a?X « x;^ + xi^ - « - <)". 

If we examine the coefficient of any power as a^'* we see that it 
is 

m^ + m^(m^ + wi, + &c.) = m^m. 

Hence the square becomes 

2m 2ma5" — 2mitw^(^/ - x^^ = 0. 

Similar expressions hold for the y and z co-ordinates. Hence 
on the whole we see that the relative kinetic energies of the several 
bodies collected at their respective centres of gravity is 

where v is the relative velocity of the centres of gravity of the 
masses m^, w,. If we express this in any kind of co-ordinates, we 
may use the Lagrangian function L to find the relative motion. 

The expression for T^ agrees with that given by Prof. Ball, but 
his demonstration is quite different. Prof. Cayley has given 
another demonstration in the same number of the Astronomical 
Notices. 

The Lagrangian function thus found may be still further 
*' modified." To avoid symbols of summation, let us consider the 
case of three particles moving in one plane under their mutual 
attractions. Let the separate masses be tn^, m„ m^, and let fi be 
their sum. Referring the system to m^ as a central mass, let the 
distances of wi„ m^ from m^ be respectively r and p, and let the 
opposite side of the triangle be R, Let the interior angle be- 
tween r and p be ^ and the exterior angle between r and R be 
X» Let be the angle r makes with some fixed straight line in 
space. We easily find 

where A = m^m^ + m^j^ + mm,JS^y 

5=:m,(mipV + ^^X)» 

O = tp- 2mmV, 

5—2 
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and V is the relative velocity of the masses m, m' calculated on 
the supposition that mj is fixed, and that the straight line r has 
no rotation round wip •Thus A, B, C are all functions of r, p, ^ 
and their differential coefficients with regard to the time. 

If we only want the changes in the form and magnitude of 
the triangle joining the three particles, we may eliminate ff by 
means of the equation 

Aff+B = Ci. 

We then find as our modified Lagrangian function 

which contains only the three co-ordinates r, p and 0]. 

24. When the geometrical equations contain differential 
coefficients of the co-ordinates f , rj, f, &c. of the system with regard 
to the time, we cannot express the co-ordinates x, y, z of any 
element of a body in terms of ^, 17, t^y &c. by means of equations of 
the form 

y =/2(f» 'ny £ &c., t), 
« ^fzit V, ?> &c., t). 

It follows, as is pointed out in our books on Rigid Dynamics, 
that Lagrange's equations cannot be employed in the form 

d dT_dT__dV 
dt d^ di" df ' 

In many of the most interesting problems in Rigid Dynamics, 
it so happens that the geometrical equations do contain ~, 

-^ , &c. For example, let a sphere be set rotating about a vertical 

diameter and be on the summit of a perfectly rough surface of 
any form. If a small disturbance be now given to it, the sphere 
may roll round and round the summit. During this motion the 
velocity of the point of contact is zero, and our mode of re- 
presenting this analytically in terms of the co-ordinates will give us 
two equations of the form 

To include such cases the equations of motion must be modi- 
fied. If L be the difference between the kinetic and potential 
energies, all the Lagrangian equations may be written in the form 

(d dL dL\ 5.^ . fd dL dL\ 5. , « ^ 
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where Sf > S17, &c. are any small arbitrary displacements consistent 
with the geometrical equations. But if these geometrical equa- 
tions be given in the form 

&c. = OJ 
these arbitrary displacements must satisfy 

5.=or 



G,S^+G^S7) + .., = 0] 

&c. 



If this were not the case, the geometrical displacement of the 
body given in applying Virtual Velocities would not be such as to 
cause the unknown frictional forces, &c. to disappear. Using the 
principle of Indeterminate Multipliers, we get 

a dL dL ^ •w,' 



d dL dL ^ J-, TT 

dtH"^''^^^'^^^^^-^^ 

&c. = 
These joined to the geometrical equations 

. = 0) 

!.=0J ' 



» • 



&c. 



are sufficient to determine the unknown co-ordinates f, 17, &c. and 
the multipliers \ fi, &c. 

It will be more convenient to write these equations in the form 
d^dJj dL ^dG dH 
dtd^' d^'^^d^^^df^"" 
ddL dL dG dK ^ _^ 

&c. = 0, 
the geometrical equations being 

= 0, H:= 0, &c. 

It is of course obvious that these indeterminate coefficients 
\, fi, &c. are merely the frictions or other resistances introduced 
into the equations in a convenient form. 

25. In order to apply these equations to the oscillations of a 
system about a state of steady motion, it will be convenient to 
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change the co-ordinates ^, ff, &a into others 0, (f>, &c. which vanish in 
the steady motion. ' Let L be thus expanded in powers of 0, <f>, &c. 
as ^explained in Art. 1, and let P and Q have the meaning given 
to them in Art. 5. 

Let us then put 

r = a+^^ V = /8 + f, &c. 

where a, jS, \, /jl^, &c. are the values of ^, nf, \,^ &c. in steady 
motion. The geometrical equations will then take the form 



&c. 



^ = 0) 
and the equations connecting 80^ B<t>, &c. will be 



O,S0 + G^S(f> + &c. = 

&c 



.=0) 

i.=or 



In these equations G^^, G^,, &c. are functions of 0, <l>, &c. Let 
a square bracket indicate that the value of the inscribed quantity 
in steady motion is to be taken. Thus [(?J means the value of G^ 
when 0, <f>, &a have all been put zero. 

The equations of steady motion may then, exactly as in Art. 3, 
be written 

&c. = 

From these the relations which exist between a, ff, &c., Xo» A*o> ^^• 
may be found. 

The equations of the oscillatory motion may be written 

with similar equations for if), '^, &c. 

26. The final determinant written for two variables 0, <l>, and 
two geometrical equations G and B in the notation of Art. 19, 
will be 
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B,y - A,, + E^ 

^'^'\j[dy'^'>ldjo\ 


B,y-A^ + E^ 


-dG' 'dH 
Idffj Idd'j 

-dor -dH 

\d<\> \ L^^ J 





^ -dG; ^ 'dH; 
'^'""Iddy'^'lde] 


B^m* -A„ + E„ 


[dOl ^ [dGl 

ld0_ "'"Lrf^'J"* 


rdGA rdGi 

cl(f> d<f> 


'dH- ^ -dB- 

Iddj^ldd']"' 


-dH' -dH' 

dO J \_d<t^ J 



0. 



It will be noticed how very much this determinant is simplified 
if the values- of \, fi in steady motion are zero. 

27. Let us apply these equations to the solution of the follow- 
ing problem. 

A heavy sphere rotating about a vertical diameter rests in equi- 
lihrium on the summit of a perfectly rough surface and being 
slightly disturbed makes small oscillations, find the periods. 

As the sphere moves about, its centre always lies on a surface 
which may be called parallel to the given surface. Let the high- 
est point of this surface be taken as the origin and let the axes of 
X and y be the tangents to its lines of curvature at 0, so that the 
equation to the suriace in the neighbourhood of is 



\Pi PJ 



Let P be the centre of the sphere, PC that diameter which is 
vertical when the sphere is in equilibrium on the summit. Let 
PJ., PB be two other diameters forming with PG a system of 
rectangular axes fixed in the sphere. Let the inclination of PC 
to the axis of Z, which is vertical, be 6, and let the vertical 'plane 
through PG make with the plane xz an angle '>^, and with the 
plane GPA an angle 0. The vis viva 22^ of the sphere will then 
be given by 

r« i (^'» + y'') + i A:« {(f + 1' cos Bf + ff' + sin'^Vr"-}. 
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Let sin d cos -^ = f , sin 5 sin -^ = 17, then we have to the 
necessary degree of approximation 

e'» + sin»^^'«-r + i7'»r 

Also let ^ + '^ = %. These transformations of co-ordinates are 
all permissible, because they do not involve any differential coeffi- 
cients with regard to the time. We thus find if i be the differ- 
ence between the kinetic and potential energies 

i=i(«'-+y-)+f{x'*-x(^»'-';r)+r+'j'*}+|(^+^. 

If a)a,, ©y, tt), are the angular velocities of the sphere about 
parallels to the axes, the geometrical conditions are 






where a is the radius of the sphere. These equations by well- 
known rulea reduce to 



X 



|r ^' sin -^ sin ^ 4- ^ coa-^ — ('^' + ^' cos d) -^ = 



y 






-^ -(^'cos-^sin^H-^ sin-ilr + (^'+ 6'cos^) - = 
a p^ 

expressing these in terms of our new co-ordinates we have 

r% 



a 



ri 



a 



The position of the system has now been expressed in terms 
of such co-ordinates, that the coefficients in the governing func- 
tions L, O, H are all constant. See Art. 2. 

The steady motion is given by x, y, |, 17 all zero, and ;^' = n. 
To find \^ fiQ we may use the equations of steady motion 

-[f]^^[i]-.[f]-°' 

where q stands for any one of the co-ordinates. Taking q = x and 
jr = y, we see that X© *== ^ ^^^ /*o ~ ^^ 

To find the oscillation we may substitute in the determinant 
and thus form the equation which gives the periods. As most of 
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the constituents of the determinants are zero, it will be more con- 
venient to form each equation directly from the standard formula 

where q stands for any one of the co-ordinates. Taking q in turn 
to be X, y, x, f j »; we find 

**(r+xV)+\=o 

*' in" - xv) + /*! = 

Putting x'=™ these with the two geometrical equations are all 
linear and ready for elimination. 

Eliminating \, /*, we have 

Substituting for f, 7f from the two geometrical equations, we 
have 



a 



► • 



To solve these put x = Xcos (jp^ + j) y = I^sin {^pt + g) so that ' 
'p is the quantity required. We obviously have 

which is a quadratic to find Jp^ 

If /3j, p^ have opposite signs the roots cannot be real, and the 
steady motion must be unstable. If p^, p, are both positive, so 
that the sphere is on the summit, the motion is stable only if 
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CHAPTER V. 



Certain suhsidia/iry determinants a/re formed from, the dynamical c?e- 
termina/nt, a/nd it is shoum that there must be at least as many 
roots indicating stability as there are va/riations of sign lost in^ 
these subsidiary determinants, and must exceed the number lost by 
an even number. Arts, 1 — 5, and 9, 

l^his is equivalent to a maasimwm and rmnimmm criterion of stability 
with similar limitations. Arts* 6 — 8, ^ 

Effect of equal roots on this test of the stability of the system. Arts, 
10, 11. 

Example, Art, 12. 

1. In order to test the nature of the roots of the determi- 
nantal equation, let us apply a method analogous to that by which 
Dr Salmon in hig Higher Algebra proves the reality of the roots 
of the equation which occurs in the determination of the secular 
inequalities of the planets. 

2. Let A be the determinant which forms the left-hand side 
of the fundamental equation, let A^ be the determinant formed 
by omitting the j)^ row and q^ column. Let A^ be the determi- 
nant formed by omitting the first r rows and r columns. Thus 
Aj = Ajj. We then have by a known theorem in determinants 

AA =A A -A A . 

It has already been noticed that if we change m into — m, the 
determinant A is changed into another determinant whose columns 
and rows are the same as the rows and columns of the first de- 
terminant. It easily follows that the minor — A^^ is changed into 
the minor — A^^ by changipg m into — m. Hence if 

Aj2 == ^ {m^) + m^ (m^) 
then Aj^ = ^ (w') — m^/r (m^) 
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Hence the product Aj^Aj^ is necessarily positive for all negative 
values of m*. 

It also follows that if A j, vanishes for any negative value of m* 
then Ajj also vanishes for the same value of m*. 

3. When the determinant A^ vanishes, we have 

AA =-A A . 

SO that A and A, must have opposite signs, or one of them must 
be zero. Consider then the series of determinants 

A A A . A 

each one being formed from the preceding by erasing the first 
row and the first column. We thus have a series of functions 
oT rrf v f hoHc ^iagr^fis regularly diminish from the n^ to the firaL 
As we may suppose the determinant A to have a row and a column 
of zeros added on at the bottom and right-hand side, but with any 
positive constant in the right-hand bottom comer, we may add 
to this series of determinants any positive constant. We have just 
proved that if any determinant of this series vanish for a n^a- 
tive value of m^ the two determinants on each side have opposite 
signs. The case in which two successive determinants vanish for 
the same value of m^ will be considered afterwards. 

We may then use these determinants in a manner somewhat 
similar to that in which we use Sturm's functions, provided no two 
successive functions vanish for the same negative value of m'. 
No variation of sign can be lost as we pass from w* = — oo to 
m" = except by the vanishing of the determinant A at the head 
of the series. And when a variation of sign is lost, it will be 
regained again at the next root, unless a root of the determinant A^ 
separates the two roots of the determinant A. If therefore in this 
passage from -co to zero, cw many variations of sign are lost 
05 is indicated by the highest power of m*, the values of m* fowrvd 
from the determinantal eqiuition must he all real and negative. 
It will also follow that the roots of each of the series of determi- 
nants are all real and negative, and that the roots of each separate 
the roots of the determinant next above it. If, however, the 
proper number of variations of signs be not lost in the passage 
from m' = — 00 to m' = 0, it does not follow that the values of m* 
are not real and negative. 

4. If the proper number of variations of sign has not been 
lost in this passage from m^^-'oo to m' = 0, this proposition 
does not leave us without information as to the nature of the roots. 
We infer that the nurnher of real negative values of m' is equal 
to or exceeds the number of variations of sign lost by an even 
number arhd unless the number of variations be even the system 
is unstable. 
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If we do not object to the labour of expanding the determi- 
nants, we might extend this theorem to determine the positions 
of the negative values of m' as well as their number. The number 
of real negative values of m* between m' = — a and w' = — /S is 
equal to, or exceeds by an even number, the number of variations 
of sign lost in the series of determinants. In this form the 
theorem resembles Fourier's theorem in the Theory of Equations. 

5. The converse of this proposition has not been proved. 
If in the passage from m" = — oo to m* = 0, the number of varia- 
tions of signs is unaltered, it is not true that the values of m' 
cannot be real and negative. Thus in the simple case 






= 0, 



where a,,, o.^a„ are aU positive, no variations are lost, yet if 

^w > Va,^ 4- Va„ the values of m* are real and negative. And if 

^18 = ^'^+^^ *^® roots are equal and negative. It will be 
noticed that the minors in this last case are not zero. 

6. In order to discover the meaning of these losses or gains 
of changes of sign, it will be convenient to make such changes of 
the co-ordinates as will simplify the dynamical determinant as 
much as possible. Let us write 

If we now change our co-ordinates by writing for 0, <f), &c 
linear expressions of some new co-ordinates, we know that we can 
clear this expression of all the terms containing the products. 
We know also that this can be done in an infinite number of 
ways. We thus have 

where the symbols 0^, 4>it &c. represent the new co-ordinates. 
Again, let us consider the expression 



the coefficients are here the values of P, ft &c. in the general 
expression for T, Art 1, Chap. IV. when 0, <f>, &c. are all put zero. 
But since T is necessarily positive for all values off, rf...^, tj,... 
it follows that TF^ is positive for all values of ff, ^'... Hence 
bjr a well-known theorem, we may by a real linear transform- 
ation of the variables clear the expression W^ also of the terms 
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containing the products ^', if>, &c. and can make the coefficients 
of the squares any positive constants we may please. We thus 
have 



2 ' 2 



It may be shown that we cannot in general clear the expres- 
sion for Tot all the terms containing 0<f>', §'<(), &c. unless C^^—C^=0, 
&c, by substituting for 0, 0, &c., any linear functions of other 
variables. As this is only a negative result of which no further 
use will here be made, it is unnecessary to supply the demon- 
stration. 

7. Using these simplifications the determinant A will now 
take the simpler form, 



A=: 



•••• > •••• > > 



= 0, 



where «» = (7^ - Ci^, &c. 

If we form from A the series of subsidiary determinants 

A, Ai,A,..., 

terminating with any positive constant, we see that when m' = — oo, 
these subsidiary determinants are alternately positive and negative, 
and when m" = 0, they become 

which are all positive if (%, a^... are all positive. Hence if 
^u> ^8a««« ^-re all positive, the proper number of changes of sign has 
been lost, and therefore the roots of the dynamical equation A = 
satisfy the condition of stability. 

If Oil, a^j, &c. are all positive, we see that V^ is a minimum for 
all variations of 0^, <^j, &c. and therefore for all variations of the 
original co-ordinates. If a^^, a„, &c. are not all positive, there will 
be as inany variations of sign lost as there are positive quantities 
in the series a.^, a^^ &c. In this case F^ is a minimum for some 
variations of 0, ^, &c. and not for others. It is shown in the 
appendix to Williamson's Differential CdkvluSy that n independent 
conditions are necessary that a quadratic expression of n variables 
should be always positive. These are given in the form of deter- 
minants and may be briefly summed up, in the statement that 
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is alwajR positive if the roots of 








A^^ + X, -a^j, 


-^I8> 


&c. 




-^12> -^M + ^> 


^.> 


&c. 




-^l«» ^88' 


^88 + ^» 


&c. 




&C. &C. 


&c. 


&c. 



= 0. 



are all real and negative. 

8. We may now put the proposition of Art 3 into another 

form. Let L be the general expression for the excess of the kinetic 

energy over, the potential energy of a dynamical system in terms of 

its n co-ordinates x, y, &c. Let this system he moving in steady 

dx dv 
motion mth constant values ^/ jt > jf > ^^* Then ifLhea maxi- 

dx dy 
mum for all variations qfx, y, &c. keeping 3r > ^x > &c. unchanged, 

then that steady motion is stable. 

If however only r of the n conditions necessary to ma^ke L a 
maximum be satisfied, then the number of roots of the dynamical 
equation which satisfy the conditions of stability is equal to r or 
exceeds r by an even number. There cannot be stability unless n — r 
is an even number. 

If the system he oscillating ahout a position of equilibrium, 

dx dti 

-jiijy &c. are all zero, and this leads at once to the condition, 

at at 

that the equilibrium is stable if the potential energy is a minimum. 

9. In this reasoning, we have for convenience excepted the 
case in which two successive determinants in the series A, A^, A,... 
vanish for the same value of w*. But this exception is of no real 
importance, for we may change these determinants into others 
whose constituents are veiy sUghtly different from those of the 
given determinants but which are such that no successive two of 
the series have a common root. In the limit, therefore, when 
these arbitrary changes of the constituents are indefinitely small, 
the roots of the series of determinants will still be real under the 
same circumstances as before, and the roots of each will separate, 
or coincide with, the roots of the next above it in the series. 

To show that these changes are possible, let us consider the 
row of determinants beginning at the last. The determinant A^ is 
a positive constant, the next A^^ is nf + a^^. Proceeding thus, 
suppose we arrive at two determinants which we may call A^ and A 
which have a common root. If we now change the constituents 
Oji, a„, ttjg, &c. into a^ + ha^, a^^ + ha^^, &c. we do not alter A^, but, 
except for the root w*= 0, we rfo alter A in an arbitrary manner. 
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When for example a^^ is altered, we alter a constituent both in the 
first row and in the first column. Since 

AA' = A A.. -A A. , 

where A' is the determinant formed from A by omitting the first 
and third rows and columns, we see that when A and A^ both 
vanish, the product AuAgi and therefore by Art. 2 both A^ and 
Aja must vanish. The determinant A is therefore altered by 
— Ag (Sajg)' which does not vanish, since A, is by hypothesis finite 
for the particular value of w* under consideration. 

If any determinant of the series vanishes when m' = 0, it is 
clear that one of the quantities a^, Ogg... must be zero. If we 
replace this by any small positive quantity, the argument will 
apply as before. 

10. It is important to consider the effect of equal roots on the 
test of stability given in Art. 8 of this chapter. 

In this case we know that the roots of the minor A^ separate 
the roots of A. If therefore A have two negative equal roots, it is 
clear that Ay must have one of them. In the same way Ajj, A33, 
&c. must also all vanish for this value of m. Since 

AA' = A A -A A . 

it follows as in Art. 2 if A and A^ both vanish, that A^^ and A^^ 
also vanish. Hence all the first minors of the determinant A 
vanish. This is the case considered in Art. 5 of Chap. I. The 
equal roots instead of introducing into the expressions for the co- 
ordinates terms which contain ^ as a factor merely render two of 
the coefficients, instead of one, indeterminate. 

In the same way if the equation Aa*0 is satisfied by three 
values of w* equal to the same negative quantity, the equation 
A^ = must have two of them, and its principal minor must have 
one of them. Reasoning as before we see that all the second 
minors of A must be zero. This is again the test that there Bhould 
be no terms which contain ^ as a £stctor. 

The presence therefore of equal roots does not in the theorem 
of Art. 8 affect the stability of the motion. 

When a system is disturbed from a position of equilibrium 
whether stable or unstable, the roots of the fundamental determi- 
nant are separated by its minors in the manner described in 
Art. 8 of this Chapter* By what has ju&t been proved, we see 
that if the fundamental determinant have equal roots, whether 
positive or negative, these do not introduce into the integrals 
terms which contain ^ as a factor. 
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11. The proposition in the last article may he made more 
general. If the mndamental determinant be reduced to the form 
indicated in Art. 7, we shall show that if A vanish for two equal 
negative values of m* which are numerically greater than the 
greatest negative quantity in the series a^^, a^, a„, &c., then 
these equal roots will not introduce any terms mto the solution 
with ^ as a factor. If a^^, a^, &c. are all positive, this reduces to 
the proposition proved in the last article. 

Following the same notation as before, we have 

AA' = A A -A^ . 

If neither A,^ nor A,, are zero, they must have the same ^ign 
when A vanishes for a negative value of m*. For their product 
is equal to A^,A„ which has been proved in Art. 2 to be positive. 
Hence all the leading first minors, viz. A,. , A„, &c. must have the 
same sign for any negative value of m^ which makes A vanish. 

By differentiation we have 

/7A 

^ = 2mA„ + o„A„ + ... - a,^„ + 2mA„ + ... - &c. 

But we have also 

A = («»' + aj Ajj + aj,m Aj, + ... 
A = -a„wA,j+ (w' + ajj A„ + ... 
&c. = &c. 

Hence if n be the highest power of m occurring in A, we have 
»» ^ = «^ + ("»•- aj ^.t + (•»*- a«) ^M - &c. 

JA 

If then A and -y— both vanish for any negative value of m* 

greater than the greatest negative quantity in the series a^^, a^j, 
&c., we have the sum of a number of quantities all of the same 
sign equal to zero. This requires that each should be zero. We 
have therefore Ajj«0, A^=0, &c. The rest of the proof is the 
same as before. 

/7A 
By differentiating the expression for -r— and substituting for 

-7-" , -j^ , &c. their values in terms of their first leading minors, 

we may extend this proposition to the case in which the funda- 
mental determinant has three equal roots, and so on. 
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12. A sphere is suspennded hy a string OA from a fixed point 
O, and is set rotating about a vertical diameter which is in the 
same straight line cw the string with an angular velocity n. A 
small disturbance is given, determins if tike steady motion is 
stahle. 

Let be the origin, and let the axis of z be vertically down- 
wards, let Ixyly^ I be the co-ordinates of -4, the point at which the 
string is attached. Let C be the centre, and let af, ariy a be the 
co-ordinates of G relative to A, Then, exactly as in Chap. iv. Art. 
22, the Lagrangian function may be shown to be 

the mass being taken unity. 

Putting x' = w> f» Vi ^'> y' fl-U zero, we see that Z is a maxi- 
mum when Xy y, f, tj are zero, the steady motion is therefore stable 
for all values oin. 

If we put fc* = I a', and m^^^ — X', so that a?, y, &c. are all repre- 
sented by terms of the form %A cos (\^ -h a), we may, by the 
methodjs of the last chapter, prove 



(^•-f)(vt»x-«).|x.. 



This equation, whatever may be the sign of ;n, has two positive 
and two negative roots. All four give stable otscillations. . 



R.A. 
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If ili/6 energy of the system be a maacimum or minimum under certain 
conditions, the motion whether steady or not is stable. Arts. 1 — 3. 

When the motion is steady, it tmU he also stable if a certain function of 
the co-ordinates called Y + W is a minimum. Art, 4. 

If there he only one co-ordinate which enters into the Lagrangian function, 
except as a differential coefficient, this condition is necessary and 
sufficient. Arts, 5, 6. 

Additional conditions when there are two co-ordinateSr Art, 8. 

1. When a system is oscillating about a position of equilibrium, 
it is well known that we may determine the stability or insta- 
bility of the equilibrium by what we may call the " energy crite- 
rion." This criterion may also be sometimes used when the 
system is oscillating about a state of steady motion. 

Let JE be the sum of the kinetic and potential energies of the 
system. Then throughout any motion of the system we have 

JS=h, 

where A is a constant depending on the initial conditions. If 0, <f>, 
&c. are the co-ordinates of the system, .B is a known function of 
0, ff, j>, <l>', &c. Suppose that some of the other first integrals of the 
equations of motion are known. Let these be 

F,{0,0',4>,4>',&c.)^C,\, 
&c. ) 

the time t being absent. 

For the purposes of this proposition let us suppose 0, 0, if>, j>, 
&c. to be separate variables unconnected with each other except 
by the equations just written down. 

If ^ be an absolute maximum or an absolute minimum for 
all variations of 0, &, &c., those corresponding to the given motion 
making E constant, then that motion is stable for all displacements 
which do not alter the cpnstants (7,, C^, &c. 
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If this proposition be not evident, it may be proved by elimi- 
nating as many of the letters as possible. If 0^ 6\ &c. be the re- 
maining co-ordinates we have 

Let h be the value of E in the given motion, and let the 
system be started in some slightly different manner so that 

E=h + hh. 

If J? be a maximum along the given motion, then any change 
whatever in 6, ff, &c. decreases E. Hence 5, d\ &c. cannot deviate 
so much from their values along the given motion that the change 
in E becomes greater than hh. 

2. Let us apply this principle to a system of bodies which 
moves in steady motion with some co-ordinates 0, ^, &c. such that 
their differential coefficients 0\ (f>\ &c. are constant, and tlie re- 
maining co-ordinates f, rj, &c. themselves constant. Let us further 
suppose that the energy is a function of ff, ^', &c., but not of 
0y <l>. By Lagrange's Equations we have the integrals 

dT_ ^ dT p o 
dff" ^' d<l>'^ ^' 

It is clear that the system can describe any one of a number 
of steady motions, which we have already called parallel motions, 
and which are determined by 

ff =p, <l>=q,&c. 

f=ar, - 97 = y3, &c. 

where p, q, &c. a, fi^ &c. are constants which satisfy all Lagrange's 
equations. 

We have thus as many relations between these constants as 
there are co-ordinates f, t;, .&c. Let the system be started with 
any initial conditions we please, then the constants 0^, (7,, ... are 
given. These being known we have as many relations between 
the constants of steady motion as there are co-ordinates 0, if>, &c. 
The steady motion is therefore determined. If the energy of this 
initial motion is nearly equal to that of this steady motion, and 
if it be a maximum or minimum as explained above, then the 
system will never deviate far from its corresponding position in 
the steady motion, and this steady motion may be called stable. 

3. Example. A top is set spinning on its point on a perfectly 
rough horizontal ground, with its axis inclined to the vertical, find 
the condition of stability. 

Let be the inclination of the axis 0(7 of the top to the 
vertical OZ; yp^ the inclination of the plane ZOO to a vertical 

C— 2 
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plane fixed in Space, and (f> the inclination to a plane through OG 
fixed in the body. Let be the apex, the centre of gravity 
which lies in OG,h==OG. Let A, A, G he the principal mo- 
ments of inertia at 0, and Jf the mass of the top. We have then 

^=f (f + t'cos^«-f~(^' + sin'^i|r'»)+%Acos^. 

By Lagrange's equations we have the integrals 

^' + '^' cos ^ = w, 
Cn CO&0 + A sin'^^' = m, 

where n and m are two constants, the former representing the 
angular velocity of the top about its axis, and the latter the 
angular momentum about the vertical. If we now eliminate 
^' and '^' we find that JP is a minimum when ^ = a, if 

G^n^ > 4iMghA cos a, 
which is the result given by other methods. 

4. The theorem of Art. 2 may be put into another form. 
Let the kinetic energy be 

T= Tee -q-^ 27^ ff<f> + ... 

Then since 0, j>, &c. are absent from the coefficients, we have 
the integrals 

2w^+ r^f + -. = Ci- Tei^ - Te,r{ - &c. 
Te^ff + 2V^f + ... = 0, - T^^ - T^t/ - &c. 

&c. = &c. 

For the sake of brevity let us call the right-hand sides of these 
equations Cj — X, (7j — F, &c. Since T is a homogeneous func- 
tion of ffy <!>, &c., we have, as in Chap. iv. Art. 21, 

2"= Tiff Y + 2V»»fV + ••• 

If we substitute in the second line the values of ff, <^', &c. 
found by solving the linear equations just written down, we have 
the determinant 

0, G, + X, G,+ Z &c. -^ 
(7,-X, Tee, Te^, &c. 
G,^Z Te^, r^, &c. 



&c. 



&c. 



&c. &c. 
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where A is the discriminant of the terms in T which contain ff. A', 
&c. This determinant is unaltered bj changing the signs of X, Y, 
&c. and is a quadratic function of C^, C,, &c., X, Y, &c. Hence 
the terms C^X, (7, Y, Sac do not occur. If then we put 



w=- 
















&c. 
&c, 
&c. 
&c. 






we have 



T^BiX+B^^rj' + &c.+ W, 



where B^^, &c. are independent of C^, (7,, &c. Now T is essentially 
positive for all values of the variables, and therefore for such as 
make C^, (7^, &c. all zero. Hence the terms involving f, rj', &c, 
are' together a minimum when ^, 17', &c. are all zero. The coeffi- 
cients B^, &c. may all be treated as constants, since f ', rj, &c. are 
all small quantities. 

K F be the potential energy, we have therefore the following 
rule. If W+Y be a minimmn for all iKiriations of f , v> &c. then 
the steady, motion i& certainly stable. It. should be noticed that 
W+ F is a function only of f, 17,. &c. the co-ordinates which are 
constant in the steady motion. 

5. If the energy be- a function of one only of the co-ordinates, 
though the differential coefficients of aU the otkers enter into its 
value, this condition is suffment and necessary. 

Let f be this co-ordinate. Then by vis viva we have 

Diflferentiating we have 

This equation must be satisfied by the steady motion repre- 

d(W+ F) 
sented by f = a. The second term ^ ,„ must therefore 

^^ . 
vanish when f = a, so that TF+ Fis a maximum or minimum. 

To find the oscillation let us put f = a + a?, we find 



Bii\ + W-\-V=h 



'2 
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where the square bracket implies that f is to be put equal to a 
after differentiation. By the same reasoning as before B^^ is 
necessarily positive, and the motion will be stable or unstable ac- 
cording as (W 4- V) w a minijnum or maximum, 

6. If we refer to Art 21 of Chap. iv. we see that this function. 
PT + F is the value of the modified Lagrangian function i' when 
^, 71 y &c. are all put zero, and the sign of the whole function 
changed. It therefore follows by Chap. V. that when TF4- K is a 
minimum the steady motion is stable. The ''energy criterion of 
stability," as far as it applies to steady motion, may therefore be 
deduced from that given in Chap. V. Art. 8. But the mode of 
demonstration adopted in that chapter gives us more information 
as the nature of the motion, while the modes of application to 
examples of the two criteria are quite different. The energy 
criterion may also be sometimes japplied to determine the stability 
of a motion which is not steady. 

[The relation between the theorem in this Chapter in which 
^= T-\- V is made a minimum to that given in Chapters iv. and v. 
in which Zf= T— V is made a maximum may be more distinctly 
perceived by the following statement. 

Let Xf y, &c. be the co-ordinates of the system, and let L be 
the Lagrangian function* so that L = T— V, then by Art. 3 of 
Chap. IV. the co-ordinates in steady motion satisfy the equations 

-i-". f-».^ w- 

Here L is expressed as a function of Xy y, &c. x\ y\ &c. 

Suppose some of the co-ordinates as 0, ^, &c. are absent from 
the expression for L, so that i is a function of ff ^ ^', &c., the re- 
maining co-ordinates, viz. ^, 9;, &c. and their differential coefficients. 
Then if we form the modified Lagrangian function as in Art. 21 
of Chap. IV, the equations (1) of steady motion become 

f.O. f = 0,te (2). 

Here, as in the Hamiltonian equations (7^, (7^, &c. ai*e the 
0, 0, &c. components of momentum, and L' is expressed as a 
function of (7^, Cg, &c. f, % &c. I', 17', &c. 

dT 
* Let u, V, &c. be the x, y, &c. components of momentum, so that u=3-,, &c. 

ax 

and let H be the Hamiltonian function. Then ff=T+ V and we easily deduce 

from the Hamiltonian equations that, in steady motion, 

dl{ ^ dll ^ , 
-7- = 0, . =0, &c. 
dx dy 

Here H is expressed as a function of », y, &c. ?/, v, kc. 
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But i'»T-C/-(7,f-&c-F 

by Euler's theorem on homogeneous functions. In steady motion 
f , 17', &c. all vanish, hence the equations (2) become 

_W^.„, -«^ = 0.4o. (3). 

where W is the value of T when f , 17', &c. are all put equal to 
zero. Here Wis expressed as a function of (7^, (7„ &a f, 1/, &c. 

It is shown in Chap. v. that if the Lagrangian function ex- 
pressed as required in equations (1) be a maximum the motion 
is stable. It is shown in this Chapter that if the function W+ V 
be a minimum the motion is stable.] 

7. To find the condition of stahility when the Lagrangian 
function is a function of two only of the co-ordinates, though the 
differential coefficients of all the others enter into its value. 

Let f, 17 be these two co-ordinates, then the modified Lagran- 
gian function as explained in Art. 20 of Chap. IV. will be a func- 
tion of f, V, ?y V oniy» 

Let the steady motion be given by f = a, 17=^, with the 
corresponding values of the other co-ordinates ^, <^, &c. Then a 
and /3 are constants. Let f = a-far, i; = )9+y, and let us expand 
the modified Lagrangian function in powers of x, y. Neglecting 
the terms of the first order, as they only give the steady motion, 
let 

Also let E— Cja— Cia. Then the condition of stability is that 
the roots of the following equation should be of the form /3>/— 1. 



B^ — A^ — Em, B^rr? — A^ 



= 0. 



If ^'^B^B^^BJ 

A = A..A. 



and = A,,B^ + A„B,, - 2A^B,, 
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be the two discrimioaDfts and the other invariant, this leads at 
once to the conclusion that the motion is stable only when 

(1) A is positive, 

(2) JS;"-e is positive and > 2 VAA'. 

If -4u 2" + -^if^ + -^M 2" ^ *^ maximum when x and y are zero, 

the two conditions are obviously satisfied. 

This condition may be otherwise expressed ; if i' be the modi- 
fied Lagrangian function, a steady motion is given by 

f -«• f "-^ f.O. ,•=„ 

This motion will be stable if for the values of f, r^ thus found, 

^^^ UW d^dfi) \d^ dv" ■*■ dfi' d^ ^didiidi'dri] 
is positive and greater than 

^ (d^n d^r d^U\ I ( d^L d^L dT \ i 
^\d^ d/rf' d^dv) 'Vr dv* d^'dv) • 

8. The nature of the motion when thus reduced to depend on 
two co-ordinates may be illustrated by geometrical reasoning. Let 
the position be defined by two co-ordinates os, y which are zero 
along the steady motion. Let these be regarded as the co- 
ordinates of a point P referred to any axes. Ihen the motion of 
P exactly represents that of the system. 

Let us construct the conies 

Then, exactly as in Arts. 11 — 14 of Chap. I v., it will be found 
convenient to transform the co-ordinates by writing 

y = a^ + \q)* 

If /A be the modulus of transformation we have /it = ajy^ — a,6j. 
It is easy to see by actual substitution, if E' be the difierence of 
the coefficients oi pq and pq^ that 

E' = iiE. 
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If the transformation be from one set of oblique co-ordinates 
to another, let a>, a>' be the angles between the axes. We then 
have 

r ^ E 

sin (o' sin a> ' 

Transforming the axes to the common conjugate diameters, the 
conies become 



-Bu'^+-B-'^=J 



r * 



the signs of a and h being such as to make these conies real. The 
equation to find m becomes 

(B^W - A^) {B^m' - A J) + E'W = 0. 

It is therefore necessary for stability that the conic a should be an 
ellipse as well as the conic b. It is also necessary that 

both roots having the same sign and the inequality being nu- 
merical 

Let OPy OP ; OQ, OQfhe the common conjugates of the two 
conies, this condition then becomes 

E area of conic 6 OQ OQf 



sina>* ^V^ " OP^ OF' 

When the system describes an oscillation with one period, Le. 
an harmonic oscillation, the path of the representative particle is 
easily seen to be 

(5,; V - ^/)/ + (5/m' - A^') 3» = constant. 

The harmonic paths are therefore ellipses. It also appears that 
the two ellipses which represent the two harmonic vibrations and 
the two ellipses a and 6 have, all four, a common set of conjugate 
diameters. 
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Any 8maM term of a high order, if its period is nearly the same as that of 
an oscilkUion of the system, may prodtice important effects on the 
m^agnittuie 0/ the oscillation. Art, 1, 

Origin qfst^h terms, with am example. Arts, 2 — 3. 

Supposing the roots of the determinantal equation to satisfy the conditions 
of stability to a first approximation, yet if a commensurable relation 
Jwld between these roots it is necessary to eocarnvne certain terms of the 
higher orders to determine whether tJiey will vltim>ately destroy the 
stability of tlie system. Art, 4. 

If a certain relation hold am^mg the coefficients of these terms, they will 
not affect the stability of the system, but only slightly alter the periods 
of oscillation. Arts, 5 — 7. 

Examples, the fi/rst taken from Lagrang^s method of fmding the oscUla* 
tions about a position of equilibrium. Arts, 8 — 9. • 

If the coefficients of the equations of motion should not be strictly con- 
stcmt, but only nearly so, the stability tvill not be affected, unless the 
reciprocals of their periods have commensurable relations with t/ie 
reciprocals of the periods of oscillation of the system. Art. 10. 

1. If we understand that a motion is called stable when toy 
small disturbance does not cause the system to deviate far from its 
undisturbed motion, it is clear that we cannot be certain of the 
stability without examining the terms of the second order. It is 
possible that some of these may have their periods so timed that 
their effects may accumulate until the motion is changed. 

Returning to the equations of Art. 3, Chap. I. we shall have on 
the right hand, instead of zero, a series of small terms of orders 
higher than the first. To find a second approximation, we substi- 
tute the values of x, y, &c. given at the end of Art. 3, in these 

terms. 

« 

They will therefore take the form iVe*** and will produce in 
or, y, &c. terms of the form ttt-t e***, where N' is of the same order 
at least as the term considered, and/ (771) has the same meaning as 
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in Chap. I. These will have to be expressed in trigonometrical 
real forms, but it is unnecessary to exhibit the process, for we see 
at once that no small term or force (whatever it may be called) of 
a high order can affect the stability of the motion unless it makes 
f{n) very nearly or exactly equal to zero. In this case its period 
is very nearly or exactly equal to one of the periods of the motion 
given by taking terms of the first order only. 

A remarkable use of this principle was made by Captain Kater 
in his experiments on the magnitude of gravity. It was important 
to determine if the support of his pendulum was perfectly firm. 
He tells us that he had recourse to a delicate and simple instru-. 
ment the sensibility of which was so great that had the slightest 
motion tdken place in the support it must have been instantly 
detected. The instrument consists of a Bteel wire the lower part 
of which, inserted in the piece of brass which serves as its support, 
is flattened so as to form a delicate spring. On the wire a small 
Weight slides by means of which it may be made to vibrate in the 
same time as the pendulum to which it is to be applied as a test. 
When thus adjusted it is placed in the material to which the 
pendulum is attached, and should this not be perfectly firm its 
motion will be communicated to the wire, which in a little time 
will accompany the pendulum in its vibrations. This ingenious 
contrivance appeared fully adequate to the purpose for which it 
was employed, and afforded a satisfactory proof of the stability of 
the point of suspension. See Phil, Trans,, 1818. 

2. Since the term Ne"^ is obtained by compounding the 
different terms in the values of a?, y, &c. it is clear that 



• • • 



n =pm^ + qm^ 

where ^, j, &c. are positive integers whose sum is the order of the 
term. It is therefore only when the roots of the dynamical equa- 
tion /(m) = are such that a linear relation of the form 

pm^ + qm^ + . . . = m^ very nearly 

exists between them, that we may expect to find important 
terms among the higher orders. The order of the terms to be 
examined will be ^ + g' + . . ., and unless this be also small, the terms 
will probably remain insignificant. If the root m^ should occur 
twice in / (m) = it is clear that the divisor / (n) will be a small 
quantity of the second order, and the term may be said (as ia 
the Lunar Theory) to rise two orders. 

3. To take an example, let us suppose a particle to be describ- 
ing an ellipse about a fixed centre of force in one focus. If dis- 
turbed it will describe a slightly different ellipse. Since 

r = a{l — ecos {nt-\- € — w) + ...}, 

=-nt-\- € -\- 2e sin (n^ + € — ft)) + . . . 
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we see that a slight change in the elements will cause variations 

in r and 6 of the period — , au additional variation in of the 

form tSn + Se and an additional variation in r of the form ha. All 
these variations should by Art, 3 of Chap. I, be indicated by 
expressions of the form 

where the values of m are the roots of the equation f{rn) = 0, 
The roots therefore of the equation /(wi) = for 8 r are m = 

and ± nV — 1, and for 8^ are w = 0, and + » V — 1. We there- 
' fore infer that any small disturbing causes of the second order 
whose periods are nearly equal to that of the particle, will cause 
important inequalities in both Sr and hd, and (since /(m) =Ohas 
two roots equal to zero) any term of long period will rise two 
orders in ZO. 

4. If the roots of the subsidiary equation are such that the 
relation 

pm^ + qm^ + . . . = w, 

holds accurately, the solution changes itscharacter. We have 
now in the value of a; a term of the form ^^3 Unless the real 
part of m^. is negative, this indicates that the system will depart 
widely from the motion which we took as a first approximation. 
We must therefore modify our first approximation (as in the Lunar 
Theory) by including in it the terms which produced these im- 
portant effects. We may then* enquire how far this modified first 
approximation indicates that the motion is stable or unstable. 
When these terms are included the equations to be solved are in 
general no longer linear, and it is sometimes impossible to find a 
solution sufficiently accurate to serve as a first approximation 
throughout the whole motion. 

5. In some cases, however, the oscillations may still be 
represented by expressions of the form 

4&C., 

where the values of n^, n^,. differ but slightly from the roots of 
the equation /(w) = 0'. Let us investigate the condition that this 
should be true, and also determine whether the changes in the 
values of m^, m,... are sufficient to affect the stability. 

Suppose that we have completed our first approximation, and 
find on proceeding to a higher approximation that the terms 

J\r e-i« + JV,e">* + . . . 
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present themselves on the right-hand side of the first set of 
equations in Art. 3, Chap. I. These terms are supposed to have 
arisen from several relations of the form 

m 

pm^ + jm, + ... = ?Wj. 

If these terms can be included in the first approximation by 
writing Wj, rij, &c. for m^, m^, &c. we have, by substitution in the 
differential equations, certain equations connecting n, M, M\ &c., 
whose left-hand sides are the same as those used in Art. 3 with 
«i written for Wj, but on the right-hand sides we have instead of 
zero the* quantities N^, N^', &c. The test of the success of the 
process is that these modifications in the values of m must satisfy 
the same relations as before. 

Now JVp if/, &c. are all at least of the second order of small 
quantities, hence up to terms of the first order the ratios M^, Jf/, 
&c. will be the same as before, so that we may put 

Mi=^Lja^y M^ — Lfi^y &e., 

following the same notation as in Art. 3, Chap. I. We also have 

Let Wj = Wi + Zm^ we find 

Similarly 

&c. = &c. 

It is evident, by the theorem of determinants alluded to in 
Art. 3, that these are symmetrical expressions. 

6. We may conveniently express these results in the form of a 
rule. 

Suppose we have to a first approximation 

Eliminate from the differential equations all the variables except 
X in the usual manner. This may be done by performing"on the 
several equations the operations represented by the minors a, a , 

&c., -7- being written for wi. Let the equation thus found be 



/(|^)aj = P,e'"^' + P/.v... 



K^^^*- ^1 
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Then all these terms caai be included in the first approximation, 
provided 

satisfy the relations 

pSm^ + qS?!!^ + ... = Sm^, 

&c. = &c. 
which exist among the roots of the dynamical equation. 

7. The general results we have arrived at may be summed 
up as follows. Though some of the terms of the higher orders 
may affect the magnitude of the oscillation, yet no term will arise 
to affect the stability of the motion unless there be some relations 
between the roots of the dynamical equation of the form 

pm^ + qm^+ ... =m^, 

where p, q, &c. are all integers. Even if such relations occur, the 
lowest order of the term is p + g+ ..., and if this be considerable 
the term will not produce any important effects until a con- 
siderable time has elapsed. If a certain relation, just found, hold 
among these terms, their only effect is slightly to modify the 
periods of oscillation, without altering the type of motion. 

8. As an exampUy let us consider a system of bodies to be 
oscillating about a position of equilibrium. We know by Lagrange's 
general solution, that the equation fijn) = is of ati even order. 
Its roots are of the form 

Whatever the numerical values of these may be we have 
^1 + ^8 + ^4 = ^1' wi2 + mg + m^ = wijj, &c. 

so that the small terms of the third, fifth, &c. orders might affect 
the stability of the oscillation. But we shall now show that they 
only affect the periods of oscillation, and not the stability of the 
system. 

Since both sides of Lagrange's equations must be of — 2 di- 
mensions in time and the impressed forces are also of — 2 dimen- 
sions, it is clear that these terms must consist of powers of x, t/, &c. 

-j^, -—, &c. and products of an even number of factors of 

ujs du 

-j2, -7^, &c. We know also, by Lagrange's solution, that the 

co-ordinates take the form 

ic = i¥j cos (at + \) + &c. 
y = J// cos (xt + \) + &c. 



f 
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The minors a, a\ &c. are also all even functions of -v- , hence 

cU 

the equation found after elimination is of the form 

/ (^j a? = P cos (at + X) + ^cos (fit-^fi) + ... 

Replacing aV— 1, — aV^, &c. by m^, m^, &c. we find by- 
Art. 6, 

P P 



Since jr(m) is of odd dimensions, and ?Wj=— w, we clearly have 
Smj=— Sm^ and therefore the test is satisfied. 

9. [As another example let us apply the rule of Art. 6 to 
some very simple case which will involve no algebraical substi- 
tutions of any length. 

The motion of a simple pendulum under the action of gravity 
may be made to depend on the equation 

• ^ + <^w=fi^ (1). 

where a and 13 are two constants and x is the inclination of the 
pendulum to the vertical which is supposed to be small The first 
approximation to the motion is 

X = M^d^^* + M^d^* (2), 

where m^ and m^ are the roots of the equation m' -f a' == 0. Our 
object is to ascertain by help of the rule given in Art. 6 whether 
the small force represented by Ba^ renders this first approximation 
unstable or merely slightly alters the numerical values of m^ and 



wij. 



The two roots are connected by the relation 

mj + wi, = (3). 

Substituting the value of x on the right-hand side, we have 

(S + *") ^ = ^^ ^i^« {M^e^i^ + M^e^) + &c. 
Hence by the rule in Art. 6 

These clearly satisfy the relation 

Sm, + Sot, = 0, 
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and therefore the first approximation taken above^ so far as the 
disturbing force Ba? is concerned, is stable. 

If the small force had been /8 (-57 1 instead of ^a?, it is easy 
to see in the same way that 

so that the relation S7i\ + Sn^ = would not have been satisfied. 
The first approximation taken above is therefore not sufficiently 
accurate to serve as a first approximation throughout the motion. 

In this example we have considered the effect of a small force 
of the third order in disturbing the stability of the motion given 
by equation (1). The same equation will obviously occur in many 
other cases of motion. For example, let a particle describe a cir- 
cular orbit about a centre of force situated in the centre. If 
slightly disturbed the equation giving the disturbance x in the 
radius vector takes the form 

•^ +a^x^px'' + yar + ... 

where a, fi and 7 are constants. Similar remarks will therefore 
apply to this case also.] 

10. When the coefficients of the equation of motion are not 
strictly constant, but yet do not vary much, then we may transpose 
the small variable parts of these terms to the right-hand side of 
the equation, and treat their products by the differential coeffi- 
cients of the co-ordinates as small quantities of the second order. 
Suppose the variable part of one of these coefficients to be p sin nt, 
where p is small, and let/(m) = be the equation givmg the 
periods of oscillation of the system when the coefficients are taken 
constant. Then it is clear that unless n is nearly equal to the 
sum or difference of two values of m, this term cannot rise into 
importance. On proceeding to higher orders we see that these 
terms cannot produce important effects unless some commen- 
surable relation between m and the roots of the equation /(m) = 
should be very nearly satisfied. 

11. It should be remarked that when the coefficients are not 
constant it is not a sufficient test of stability that they should 
always satisfy the conditions of stability obtained by giving them 
their instantaneous constant values. Thus if the equation of 
motion were 

rf^ 1 



^ + i;2^ — ^' 
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the coeflScIent of x is always positive, yet as the equation is satis- 
fied hj x^ahJtyX may become as great as we please. 

Even if the coefGcients are nearly constant, we must yet ex- 
amine, by the rules just given, if their small changes are so timed 
as gradually to increase the oscillation until the divergeiice from 
the given motion is no longer small. 

[Suppose a system to be oscillating so that its motion is de- 
termined by the equation 

where y is a known function of t, which during the time under 
consideration always lies between ^ and )8"^ the latter being the 
greater. Let the system be started with an initial co-ordinate x^ 
and an initial velocity x^ m a direction tending to increase x. 
It may be shown that the system will begin to return, i.e. x will 

begin to decrease before x becomes as great as w a?/ + -^ , If 
+ m, Twi' be two successive maximum values of x, we may also 
show that rri cannot be so great as ^ m, and that the time from 

one maximum to the next lies between 7: and 7^ 1 



H. A, 



CHAPTER VIII. 



The ffamiltonian Ckaracteristio or Principal functions when/ound deter- 
\ mine at once the motion of the system from one given position to 

arwiher, and whether the motion is stable or unstable. Arts, 1 — 3. 

Examples with a mode of effecting the integration S= JLdt in small oscil- 
lations. Arts. 4 — 7, 

The Characteristic function supplies the condition that tfte motion is stable 
cw to space only^ while the Principal function gives the conditions that 
it is stable both as to space and time. Art. 8. 

In wJiat sense the motion is unstable if either of the two ffamiltonian 
functions is a minimum. Arts. 9 — 14. 

1. If we had any convenient methods of finding the Hamil- 
tonian Characteristic or Principal function, we might determine 
without diflSculty the conditions of stability of a dynamical system 
at the same time that we deduce the integrals of the equations of 
motion. But it is very difficult to discover either of these func- 
tions by an d priori method. We have indeed differential equa- 
tions which they must satisfy, and Jacobi has taught us what kind 
of solution will serve our purpose. But the difficulty of finding 
these solutions is as great as that of solving the equations of 
motion. For these reasons it does not seem necessary to dwell on 
the uses of these functions. 

2. Suppose the Principal function S^ of a dynamical system to 
have been found in terms of the initial co-ordinates 6^, (f>^ and the 
co-ordinates 0, (j> and the time t. Let the semi- vis viva be given 

where P, Q, R are known functions of 0, (f>. Let P^, Q^, R^ be 
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the values of thesa when 0^, (f>Q are written for 0, <f}. The final 
integrals of the equations of motion are then given by 






#0 

Let the system receive any disturbance at the time ^ = 0, so 
that while starting from the same initial position, its initial ve- 
locities are slightly altered. Let a?/, y^ be these initial changes 
of 0^' and 0o' and let + x, <f>-\-y be the cO-ordinates of the 
system at the time t Then we have 

(PS (P8 I, . , n ' 



d0,d0 d0,d(t> 

^d<f^0'"'dm^^^^''^^'^'i 

Here x^, y^ are small arbitrary quantities, hence x and y will 

!.« « 11 -4? c xu X- r <^S d^8 d^S d'S 

be small if none of the ratios of ,-77-^ , ^^ ,. , ,. ,^ .or 



d0^d0' d0,d(l>' d<l>^d0' d<t>^d4> 
to the determinant 

d'8 d^S ^S d^8 



d0^d0 d(t>^d(j> d0,d<l> d<j>^d0 * 
be large. 

If the initial position as well as the initial motion be altered 
we may find, by a precisely similar process, the conditions that 
X and y should be small. If the system have more than two inde- 
pendent motions, we have more than two co-ordinates, but the 
conditions of stability are found in the same way, 

3. If aj and y be small throughout the whole motion from the 
one given position to the other, not only does the system not 
deviate far from its undisturbed course, but the system at any 
instant is also very nearly qoincident with its undisturbed place 
at the same time. It is important to notice this, for the word 
" stability " is sometimes used in a diflferent sense. 

This condition of stability may be put under a form in which 
no reference is made to 8, Let u and v be the components of 
momentum of the system corresponding to the co-ordinates 0, <^ 

respectively, ie, let ^ = -5^, , v— jT7> and let these be expressed 



^> 
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as functions of 0^,if>^, 6,^ and t Then the preceding equation 
may be written 

du dv 

du dv ^ 

where a and /8 are small arbitrary quantities. The condition of 
stability is that the values of x and y thus found should be small. 

4. As an example let us consider the case of a projectile. 
If ^ be the horizontal, and (f> the vertical co-ordinate of the par- 
ticle, we have 



The equations to find x and y are evidently 



1 



* * 



Hence the system continually deviates more and more from its 
undisturbed place. 

5. In order to calculate the form of 8 when a system is oscil- 
lating about a state of motion, it is convenient to choose as co- 
ordinates some small Quantities x, y which vanish in the given 
state of motion. Let the Lagrangian function be written in the 
form 

where L^ is a homogeneous function of x, y, x\ y\ Then by a 
theorem of Euler's 

where the S's imply summation for all co-ordinates. 

As in Art. 9 of Chap.. IV. we have 

d dL^ __ dL^ 
dt dx dx 

and the oscillations are given by 

d dL^ __ dL^ 
dt da/ dx 
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Hence we find 
Integrating we have 

«-/4<«..[»§]>is[.§];. 

Thus the integration has been effected, but in order to express 
S as a function of a?o, y^^ x, y and ^, it will be necessary to find 
^o'> y^y ^> y' ^ terms of these quantities. 

6, As an example, let the position of the system depend on 
one co-ordinate x and let 

where the coefficients are all constants. We then find by the 
process just indicated that -4^ = and 

m^„ {a? + x^) [e ^ + e'^) - 4ixx ^ 

A 
where m? = -^ . Applying the criterion of stability we find that 

11 

, , will finally become small if m is reaL The motion is there- 

CbJCJXiC 

fore unstable or stable according as A^, ^^ have the same or 
opposite signs. 

7. If the position of the system depend on two co-ordiniitcs 
a, y, let 

i = i, + JjSJ + ^^ + B^x' + Bjf' 

+ Mua? + ^« a^ + i ^«»* + i ^u*" + B,^x'y' + i B„y'* 
+ C^xx' + C^ + C^x' + C^y\ 

We then find 
S= L,t + B, («-a;J + B, (3^ -yj 



where 



<r = :?" {xx' - xX) + % (-y + a'V - ic,y: - a^.'y.) + ^"^ (yy'- y^.'). 
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If we now express a;', y\ a?/, yl in terms of a?, y, x^ y^ and f, we 
find for (T a fraction whose numerator is a homogeneous quadratic 
function of a?, y, x^ y^ the coefficients being Imear functions of 
exponentials of ^, and whose denominator is another linear function 
of the same exponentials. These exponentials become sines and 
cosines when the motion is stable. Thus when the given motion 
is steady the simplest inspection of the form of a will determine 
whether the motion is stable or not. 

Referring the motion to principal co-ordinates for the sake of 
brevity, and writing 2 (? = (7^, — 0,^, we find that a must satisfy 
the diflferential equation 

This equation is obviously satisfied by such a function as that 
just described. The solution of this equation may be reduced to 
linear equations and thus <t may be found. But it is unnecessary 
to dwell on this, for this would be equivalent to returning to the 
Hamiltonian equations* 

8. If we wish to determine the condition that the general 
course of a dynamical system is stable without requiring it should 
be near its undisturbed place at any the same time, it is more 
convenient to use the Characteristic function. Suppose that the 
Characteristic function has in Jacobi's manner been expressed as 
a function of the co-ordinates 6y <f>y the constant h of vis viva and 
two arbitrary constants a^, a^. Then 

V=f{0,<f>,h,a,) + a^. 

The relation between and (f), which may be called the equa- 
tion to the path of the system, is given by 

dr_ , 

where \ is another constant. Let the system be disturbed from 
the same initial position so that the whole energy is unaltered. 
The change in ^ corresponding to any given value of is found 
from 

Let A be the initial value of -^— s , then 

d^V 



-A 



da^d<^ 
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The condition that the path should be stable is that the coeffi- 
cient of Sa^ should not be large. 

We might also use the function called Q by Sir W. R Ha- 
milton, but it seems unnecessary to dwell more on this subject. 

9. The instability of a system may be deduced from the 
Hamiltonian Principal or Characteristic fiinctions, expressed as a 
minimum. Suppose a dynamical system to move from one posi- 
tion A to another J9 in a time t, then the motion may be found by 

making the first variation oi 8=^ I Ldt equal to zero, the time of 

transit being constant. The constants of integration are deter- 
mined by the conditions that the co-ordinates have given values 
when ^ = Oand t=t To determine whether /S is a maximum or 
minimum or neither we must examine the second variation and 
here we have the assistance of Jacobi's rule. The determination 
of the constants will depend on the solution of equations and may 
lead to several different kinds of motion from A to B, One of 
these will be the actual motion. Let us move B along this until 
one of the other motions coincides or as we may say approaches 
indefinitely near to this actual motion. We have then reached 
a boundary beyond which the integration must not extend if ;Si is 
to be a maximum or minimum. See Todhunter's History of the 

cPL 
Calculus of Va/riationSy page 251. Further -jmy if ^ be a co- 
ordinate, is positive throughout the limits of integration, so that 
8 will be a minimum and not a maximum. 

10. When there are several co-ordinates 0, (f>, &c. which are 
to be found as functions of the time, we may easily show that 
Jacobi's condition is a necessary one, and this is all that we require 
for the next proposition. If the system can move in two ways 
from A to B, then 88=0 along each, and therefore when these 
two ai'e adjacent we have both SS = and S{8-{-BS) = 0. This 
shows that the second variation can be made to vanish by taking 
one variation through the other. This second variation will then 
be the same as the quadratic term of the series obtained by 
changing the co-ordinates 6, <f> into ^ + S^, ^ + S<^, because we can 
take ^6 = and S*^ = 0. Hence as the sum of the terms of the 
third order does not, in general, vanish for this displacement, it is 
clear that 8 cannot be either a maximum or a minimum. 

Let the actual motion be from A to J5, and let a neighbouring 
motion starting from A lead the system to a position G reached in 
the same time along the actual motion before reaching jB. Then 
we can show that a variation of the actual motion from J. to JS 
can be found which makes 8*5 of any sign. Let P be any position 
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on the neighbouring motion before reaching O, and Q, one on the 
actual motion after passing C. Then considering P and Q as fixed 
and also the time of transit, the motion along PCQ cannot 
make S5= 0; for this condition is known to lead to the ordinary 
dynamical equations, and it is clear that (impulsive forces being 
set aside) no actual motion can be discontinuous. But there is 
discontinuity at (7, for otherwise when the system is started from 
B towards A, two courses would be open to the system on arriving 
at a Hence the first variations of 8 tor an imaginary motion along 
FCQ are not zero, and therefore may have any sign. But since 
the discontinuity at is of the first order of small quantities, this 
first variation is of the second order. Now the value of 8 for the 
actual motion is equal to that along the neighbouring motion to 
C and then along the actual motion to B. Hence, P and Q being 
still fixed, variations of the actual motion from ^ to -B can be 
found which make S^8 of any sign. 

11, Let us apply this theory to determine the stability of a 
given state of motion. First let us suppose the given motion to 
be steady and to depend on only two variables. If we use the 
function 8 there will be one co-ordinate and the time, if V two 
co-ordinates. Let the system be disturbed at any moment by an 
alteration of the velocities of its several parts, so that the initial 
position of the disturbed motion is an undisturbed position. If 
the motion be stable the system will oscillate about the un- 
disturbed motion, the oscillation repeating itself at a constant 
interval. It follows therefore by Jacobi's rule that 8 or V cannot 
be a minimum for a period longer than the time of a half-oscilla- 
tion. If therefore 5 or F be a minimum for all variations, starting 
from A and ending at B, where 5 is a position on the steady 
motion reached by the system at an ijiterval as long as we please, 
then the motion is unstable. 

If we give a meaning to the word ** stable" somewhat different* 
from its usual signification, we may extend this proposition to 
determine a test of the stability of any motion, whether steady or 
not. All we have assumed is, that, if the motion be not altogether 
unstable, there are some disturbances which will cause the system 
periodically to assume the same positions as it would have done if 
it had been undisturbed, but the interval of these periods may be 
any whatever provided the first be finite. If we use the Charactert 
istic function, these disturbances must be such as not to alter the 
constant of vis viva, and if the Principal function, they must be 
such as to bring the system to an undisturbed position in the same 
time. 

* This meaning does not always agree with the results of Art. 14, Ohap^ iv* 
[See also Arts. 17 and 18, Chap, iv.] 
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12. Next let us suppose that as the system proceeds from the 
initial position A along the actual motion, 8 ceases to be a 
minimum at some position B, The conditions for a minimum 
are of two kinds. Suppose the system to depend on two co- 
ordinates 0y (f>, and let L be the Lagrangiah function, then (1) we 

must have -wg and ;jB^ jTts "" TtZTtT? ^*t positive, and (2) it 

must be possible to choose three arbitrary constants which enter 
into a very complicated expression, so that this expression may 
never become infinite between the limits of integration. The first 
condition is clearly always satisfied since the vis viva of any 
system is necessarily positive for all values of 0' and (f>. The 
second condition will fail if there are two neighbouring motions by 
which the system can proceed from A to any position between 
A and B, If this be the mode of failure, it is clear from the 
reasoning of Art. 2 that the conditions of stability are satisfied for 
one kind of disturbance, and that therefore some at least of the 
harmonic motions are stable or oscillatory, thoitgh the motion may 
be unstable for a different kind of disturbance. 



13. [These conditions become much simpler when the position 
of the system is determined by one co-ordinate, or when the 
Lagrangian function can be reduced to depend on one co-ordinate. 
Let this co-ordinate be so chosen that it vanishes along the given 
motion, and let us also suppose that both it and its differential 
coefficient with regard to t, are small for all neighbouring con- 
strained motions. Let this co-ordinate be called 6 and let the 
Lagrangian function be 

L = L, + A,0 +S,ff+^ A,^e^ + 1 BJ^ + CJff. 

Then since the Lagrangian equation of motion is satisfied by 
hypothesis when ^ = 0, we have A^ = B^ where the accent, as 
usual, denotes differentiation with regard to t. 

If the system be now conducted from the initial position A to 
any other position By both on the given motion, by any neighbour- 
ing mode of motion, we have 

If d = u be any solution of the Lagrangian equation 

B. A. 8 
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we may write the function 8 in the form* 

The second term is essentially positive, since B^^ must be 

positive. Hence iS is a minimum along the given motion unless 

d 6 
we can so choose the arbitrary displacement d as to make -^ - = 0. 

Cut U 

This gives O^cu where c is some constant. But must vanish 
at the two limits A and By hence this choice of is excluded 
unless there is some neighbouring mode of motion by which the 
system could move freely from the given initial position A to the 
position B. The result is that 8 cannot cease to be a minimum 
before the first instant at which some neighbouring motion will 
bring the system (starting from A) into coincidence with some 
contemporaneous position on the given motion. If the given 
motion be steady, it follows that 8 cannot cease to be a minimum 
before a time which is half that of a complete oscillation. 

We thus have a test of stability. If the system depend on one 
co-ordinate and if ^ he a minimum when the limits of integration 
are from the initial position A to ail positions oh the actv>al motion, 
that motion is imstable. But if S cea^e to he a minimum at some 
point C, tlien the actual motion is stable from A to C] 

* [Following Lagrange's rule we may write the second term of 5 in the form 

The quantity outside the integral sign is to be taken between the given limits 
and is zero, since $ vanishes at each limit. Let us now put 

It is dear that this value of X cannot be infinite between the limits of integra- 
tion unless u vanishes. For by hypothesis u and u' are both finite and the co- 
efficients B^i and (7^ in the Lagra^gian function are also finite. It then easily 
follows from the equation 

that (0'u + 2X)»=Bii(Ju + 2V). 

Hence the second term of 8 becomes 



lj^n{e'-eipt. 



which is the result in the text This might also have been deduced from Jacobi's 
general transformation with one independent variable given in Prof. Jellett's 
Calculus of Variations or Prof. Price's Differential Calculus, 

If u=0 between the limits of integration this transformation fails, but, as is 
evident from the argument in the text, we choose e=uio represent a neighbouring 
free motion such that u=:0 just before the system reaches A ; also B is so placed 
that the next instant at which u=0 is after the system has passed B.] 
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If the angle A OB = ff, we therefore have 
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^ <f>{a) 

If the neighbouring path be a free path desciibed with the 
same energy, its equation is 

p = X sin p$, 
where ia measured from the radius vector OA, and X is an 
arbitrary constant This free path will cut the circle again in 
some point 0. If the angle AOO = y, we have py^ir. 

If B coincide with G, we find by substituting this value of o in 
the expression for V, that the second term of Fis zero. If 5 be 
beyond C but such that the angle GOB is less than 7, draw two 
free paths one from A aa before and the other backwards from B 
to meet the former in some point P. Then the angle AOP=^ff. 
If the particle be conducted from ^ to P along one path and 
from P to B along the other, we find that the excess of the 
action over the action in the circular arc AB is equal to 

'■ ^i'sin^^. 

Since p$ is gi-eater than tr and less than 2ir, this excess is 
negative. The action along the circular arc is therefore not a 



>e^r^s^=^s^iHB^M«59i 
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minimum if 5 be beyond the first intersection of a neighbouring 
free path. 

Lastly we may show that the action is a minimum if B lie 
between A and C, To prove this we write the integral in the 
expression for V in the form 

The first of these two terms is zero since p vanishes at each 
limit. Following Lagrange's rule we make 

The integral then becomes 

This is always positive and the action along the circular arc is 
a minimum. The argument however requires that \ should not 
become infinite between the limits of integration. It is easy to 
see that 

\ =p tan (p0 + E), 

where E is a, constant to be chosen at our pleasure. But if fi 
exceed - it is impossible to choose E so as to keep X finite be- 
tween the limits of integration. Hence the action is a minimum 
only if the angle A OB subtended by the limiting positions at the 

centre is less than - . . The circular orbit is therefore stable if ©* 

. . P 

be positive. 

If, however, p^ be negative, the expression for \ changes its 
character. Writing — g" for jp*, we find 

^^^Ee-^, 
q—X 

thus the value of \ can be chosen so as not to be infinite for aU 
positive values of 0. In this case the function F is a minimum 
for all arcs AB however distant B may be from A. The circular 
motion is therefore unstable.] 
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